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DEFLECTIONS OF STRUCTURES IN THE INELASTIC RANGE 


Kurt H. Gerstle,! J.M. ASCE 
(Proc. Paper 1290) 


SYNOPSIS 


An approach to the problem of finding deflections of structures of perfectly 
plastic material in the inelastic range is developed. Charts relating loads and 
deformations of rectangular beams are furnished, and may be used to analyze 
indeterminate structures composed of members of such cross sections in the 
inelastic range. 


INTRODUCTION 


Two basic approaches to the problem of finding deformations of structures 
in the inelastic range have been taken in the past. On one hand, certain con- 
venient approximate methods2 are available, which, due to the simplifying as- 
sumptions involved, are of questionable accuracy for certain cases. On the 
other hand, rigorous methods have been devised, which, however, are of such 
complexity as to make them unsuitable for many applications. Fh 

This paper attempts to simplify the determination of deformations and the 
analysis of moments of indeterminate structures in the inelastic range while 
subjecting itself to the same assumptions as the rigorous approach mentioned 
above. 

The integration of unit angle changes results in expressions for total angle 
change and tangential deviation which may be written as dimensionless coeffi- 
cients dependent on the loading condition, to be multiplied by a factor depend- 
ent only on the properties of the member. The coefficients are plotted for 
different loading conditions of members of rectangular cross section. Using 


Note: Discussion open until December 1, 1957. Paper 1290 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 3, July, 1957. 

1. Instr. of Civ. Eng., Univ. of Colorado, Boulder, Colo. 

2. A summary of various methods of inelastic deflection calculations is con- 
tained in Progress Report No. 9, “Plastic Strength and Deflections of Con- 
tinuous Beams,” by K. E. Knudsen, C. H. Yang, Bruce G. Johnston, and 
Lynn S. Beedle, Fritz Engineering Laboratory Reprint No. 86. 
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the graphs, the deformations of a structure due to a given moment diagram 
may be determined. 

Statically indeterminate structures may be analyzed in the inelastic range 
by assuming a statically compatible moment diagram and performing a check 
upon the conditions of geometry, utilizing the furnished graphs. 

The application of the proposed method to the analysis of structures com- 
posed of members of I-section is discussed. _ 


Assumptions 
. Idealized flat-topped stress-strain diagram. Strain hardening is neglected. 
Equal properties in tension and compression. 


2. Plane sections remain plane after bending in both elastic and inelastic 
ranges. 


Small deflections. 
. Internal shear, axial forces, and possibility of buckling are neglected. 


5. Effects of residual stresses, stress reversal, and sequence of loading are 
neglected. 


Since the neglect of strain hardening results in error on the side of safety, 
whereas the assumptions stated in items 3,4, and 5 give results on the unsafe 
side, therefore the errors will tend to compensate. 


Inelastic Bending of Rectangular Beams 
It is desired to find a relationship between the rotation of adjoining sec- 


tions of a beam, and the applied moment. In the elastic range, such a rela- 
tionship is furnished by the well known expression 


In the inelastic range, considerations of statics and geometry help to establish 
a similar relationship. Equating external and internal moments on any sec- 
tion of the rectangular beam, 


t (c7-3y?) . 


From geometric relationships, 


Cc Ey 


Solving the two conditions simultaneously, the required relation is seen to be 
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in the inelastic range. 
Angle changes and tangential deviations of flexural members may be deter- 
mined by direct integration of the unit angle changes: 


el? 


sionless coefficient C dependent only on the ratio of moment to the yield point 
moment Myp: to be multiplied by a factor involving only the properties of the 
member 


i | where @ap is the total angle change between Points A and B, and 6 ap is the ba 
a | tangential deviation of Point A from a tangent from Point B, if x is measured es 
4 from A. The results of these integrations may be written in form of a dimen- . + 
| 


L2 
YP for tangential deviations, 


and Cy EI for angle changes. 


Since the relationship between load and resulting deformation is non-linear 
in the inelastic range, the coefficients C1 and C2 must be determined sepa- 
rately not only for different types of load and points of application, but also 
for different magnitudes of inelastic loads. 


Deformation of Beams Subject to End Moments 


In this section 2re presented procedures and results for evaluating the 
coefficients Cj and C2 for members subjected only to end moments. Such 
cases occur in structures whenever an unloaded member helps carry a load 
imposed on an adjoining member; unloaded spans of continuous beams and 
legs of rigid frames may be cited as examples. Furthermore, the results of 
this section will be used in the evaluation of deformations of beams subjected 
to concentrated loads. 

The moment variation in such beams, according to Fig. 1 b, is 


(h- 28 Ky) 


therefore, according to Eq. 3, the total angle change is given by 
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and the tangential deviation, according to Eq. 4, is given by 


a,  Dellected Shape 


b, [Moment Diagram 


Fig. 1. Beam subject to End Moments. 
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The coefficients Cy and Cg are given by the results of the integrations in 
Eq. 5 and 6. Figures 2 and 3 furnish the values of these coefficients for all 
combinations of end moments hMyp and KMyp. Tangential deviation dpa is 
easily determined from the previous two quantities by the geometry apparent 


in Fig. 1 a: 
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Deformation of Beams Subject to Concentrated Loads and End Moments 


For purposes of finding deflections of beams subjected to concentrated 
loads, as well as to satisfy conditions of geometry for continuous beams and 
frames, it is necessary to determine the end slopes of such members. In this 
section a method for doing this will be set forth, making use of the coeffi- 
cients for angle change and tangential deviation found in the preceding section. 
The method presented is restricted to the case of one concentrated load on 
any one span. It could easily be expanded to include the case of any number 
of loads at any location, but it should be remembered that each combination 
requires an entirely new set of calculations, since the superposition of effects 
of individual loads is not admissible in the inelastic range. 

Since information for beam segments with straight line moment diagrams 
has been obtained in the preceding section, it is convenient to treat the two 
parts of the beam, shown in Fig. 4, AC, of length aL, and BC, of length bL, 
separately, each one as a beam with cranked-in end moments, kaMyp and 
hMyp, and kRMyp and hMyp, respectively. Given the end moments KAMyp and 


hy 


a, Deblected Shape 


Al Typ 


6b, Moment Dia gram 


Fig. 4. Beam subject to Conc. Load and End Moments. 


kBMyp, the magnitude, and the location of the load P, the moment at the point 
of load application may be found by statics: 


hl (bk, + ak, + 
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A study of the deflected curve of Beam AB in Fig. 4 a will lead to the fol- 
lowing useful relationships: 


CA 


doc +¢& 


and doe + 9. 


The end slopes are now given in terms of the tangential deviation and total 
angle change of the individual beam segments AC and BC, the values of which 
may be determined by means of the coefficients given in Figures 2 and 3. De- 
noting the coefficients pertaining to segment AC by the subscript a, and those 
pertaining to segment BC by subscript b, at the same time adhering to the 
previous subscripts 1 and 2 to denote tangential deviation, and angle change, 


respectively, there result the general formulas for the end slopes of the type 
of beam considered here: 


10. 


and - + b Cy ) Let 


Deformation of Beams Subject to Uniform Loads and End Moments 


The same principles employed in the previous sections may be used to de- 
termine the deformation of beams subjected to uniform loads and end mo- 
ments. To avoid extremely long and cumbersome formulas, the beam will be 
separated into two segments. These segments will extend from each end of 
the beam to the point of maximum moment, at which point the origin of the 
variable of integration will be taken. One segment AC of the beam as de- 
scribed, and shown in Fig. 5, will be considered; the total angle change be- 

ML 
tween A and C, 9 ac given by Co = , and the tangential deviation 


EI 
M_ L2 


given by Cy Tt » are to be found; the moments kMyp and hM,,, and the 
load intensity w are known. These three values are not independent of each 
other, if hMyp, as assumed, is the maximum moment of the Beam AB. Given 


KMyp and hMyp, the load intensity w may be found by statics, according to 
Fig. 5 a: 


2[h-k) . 11. 


The moment variation in terms of x measured as shown is then 
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a, Deflected Shape. 
My Allyp q 


k "yp 


b, Moment Diagram 


Fig. 5. Beam Segment under Uniform Load and Moments. 


From previously presented theory, the total angle change between ends of the Pi 
segment is given by Wg 


4. + al, 7 


and the tangential deviations are - 
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In these expressions the integrand has been arranged in non-dimensional 
form, so that integration results in the dimensionless coefficients C; and Co, 


M_L2 ML 
to be multiplied by the familiar factors —F— , and —F- 


Before the indicated integrations can be performed, it is necessary to lo- 
cate the limits of the elastic and inelastic portions of the beam segment. This 
may be done by setting the moment expression equal to tMyp and solving 


for x: 


All the information is now assembled to evaluate the desired angle changes 
and deviations for the beam segment. Values of coefficients Cy and C2 are 
presented in Figures 6, 7, and 8 for all possible combinations of the end mo- 


ments KMyp and hMyp.- 


ks 


, Moment Diagram 
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ec, Freebody 


Fig. 9. Beam subject to Uniform Load and End Moments. 
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Having evaluated the deformations of each of the two beam segments AC 
and BC, of the beam AB shown in Fig. 9, individually, they may be combined 
to arrive at the required end slopes of the beam. Fig. 9 a shows the relation- 
ships of the deformed shape of the beam; these will be recognized as being 
identical with those of Fig. 4 a, and are therefore given by Equations 10. 

If for any beam the end moments kaMyp and kgMyp, and the load intensity 
w are given, there remain two items to be evaluated Talore the above Equa- 
tions 10 can be solved: the maximum moment hMyp and its location. Both of 
these can be found by statics. Referring to the free body shown in Fig. 9 c, 
the general expression for the moment in terms of the variable x measured 
from End A of the beam is written: 


(1+ bake 


Differentiating with respect to x and setting the result equal to zero allows a 
solution for the distance aL of the point of maximum moment from End A: 


al [2 She] L ‘ 
bh [4 + [ky ] 1 


Resubstituting the value for aL in the general moment expression furnishes 
the value of the maximum moment hMyp: 


ky + ke lk wl? 


and therefore 


To summarize the steps taken in the evaluation of end slopes for uniformly 
loaded beams in their correct order: 


1. Calculate a and b, Eq. 16. 

2. Calculate h, Equation 17. 

3. Find the values of coefficients Cjg, Cip, C2q, C2b for beam segments AC 
and BC from Figures 7 and 8. 

4. Evaluate end slopes 9,4 and 6p, Equation 10. 


If extensive computations involving uniformly loaded beams are required, it 
L 

is convenient to plot the dimensionless coefficients of ie in Equations 10 

for direct evaluation of end slopes. A different graph is required for each 

load intensity. 


Application of the Method to the Finding of Deflections 


Graphs 2, 3, 6, 7 and 8 furnish information for the evaluation of angle 
changes and tangential deviations for portions of beams and frames subject to 
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Fig. 10 : Example Problem l. 


Simply Supported Beam under 


Concentrated Load. 


Deflected Shape and Moment 


Diagram. 


037 147 7 -2]] +223 17 O73 
50 Li .37 .37 56 .56 .04! 164 -3lo 374 +25! 19 .063 .082 
55 122 41 Al 62 62 .04% .2] 44 +26 2 09 


60 133 47 47 Jo Jo 052 209 253-30 466 +309 23 a 
65 145 57 57 82 82 .063 254 273-464 546 +385 25 .083 127 
Table 1 : for Prob.1l. 


Fig. 11 : Plot of Deformations a 
vs. Load for Example Problem l. i 
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either concentrated or uniform loads. Anyone familiar with the calculation of 
elastic deflections by means of the well known moment area method should 
have little trouble in using these charts to find any desired deformation. As 
in the moment area method, it is essential to draw a picture of the deformed 


shape of the structure, indicating the geometrical relationships needed for 
the calculation. 


Example Problem 1: 


Given: Simply supported beam under concentrated load as shown in Fig. 10. 


Required: Curves of End Slopes @,4 and 6p and Deflection under load 4c vs. 
load up to collapse. 


To find the end slopes 6,4 and 6p, Equations 10 are used, which for the par- 
ticular position of load are written as 


/ 
(95a - - Gq) , 


4 2 
and 65 -(5Ga - 9 Cr + 5 Cy 
Coefficients Cjg and Cjp are found from Fig. 2, and Coefficients Cgq and C2p 
from Fig. 3. To find the Deflection Ac, relations obtained from the deflected 
structure of Fig. 10 are used: 


2 


Cia refers to the deviation of Point C from the tangent through Point A and is 
obtained from Fig. 2. 


Application of the Method to the Inelastic Analysis of 
Indeterminate Structures 


In regard to the analysis of indeterminate structures in the inelastic range, 
these points are to be borne in mind: 


1. Conditions of statics apply under all circumstances. 
2. Conditions of geometry may be used if the deformations are obtained ina 
manner Similar to the one described in the preceding sections. 


3. Principle of superposition is invalid. The total effect of all loads must be 
lumped together. 


The procedure used here for analyzing st.uctures is to assume a statically 
compatible moment diagram, and to check whether the prescribed conditions 
of geometry are satisfied. If they are not, then the assumed moments must be 
revised till both conditions of statics and geometry are fulfilled. This pro- 
cedure is explained in detail in Example Problem 2. 


Example Problem 2: 
Given: Rigid Frame under vertical load as shown in Fig. 12. 


Required: Vertical deflection under the load Ag and side sway deflection Ac 
up to collapse. 
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The frame shown in Fig. 12 and 13 is estimated simultaneously, such that 
two conditions of geometry are satisfied; using the convention that clockwise 
rotations, or deflections connected with clockwise rotation of chord, are 
called positive, these conditions may be written, referring to the deformed 
structure shown in Fig. 13, as 


42a = 
and B or i, + 


Moments Ma and Mg are guessed at, the other moments are calculated by 
statics, and a check is performed whether the above conditions are fulfilled. 
This work is shown in Table 2. An elastic analysis was first performed to 
get an idea of the relative magnitudes of moments, and for further work it 
was recalled that the smaller moments always tend to catch up with the 
larger ones. 

A plastic hinge forms under the load when it reaches a value of approxi- 

M 


mately 10 =". After this instant, the structure becomes indeterminate to 


the first degree, and only the second of the above mentioned conditions of 


geometry must be satisfied. At this time, two requirements should be 
observed: 


1. The graphs may not be used to determine angle changes through a yield 
hinge. 


2. No line tangent to the deflected beam at a point where a yield hinge is 


located may be used, since the inclination of such a tangent is indetermi- 
nate. 


To find the value of 6¢ for the above condition of geometry without vio- 
lating these rules, the fact may be utilized that the vertical deflection of the 
plastic hinge under the load must have the same value whether measured 
from the right or from the left. Thus, 


or - [dy + dec). 


Substituting this value of 9¢ in the above condition of geometry, it will take 


the form 
-[% [eat See) he 


All quantities in this expression may be calculated without violating the stated 
two rules. 


Table 2 shows all calculations connected with this work. The required re- 
sults may be taken directly from this table. 
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Example Problem 2. 
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Application of the Method to Different Cross Section Members 


The entire preceding discussion has been concerned with beams of rectan- 
gular cross section. It remains to mention the applicability of the suggested 
method to the analysis of beams of different cross section. 

Above all, it should be noted that the basic method may be used for mem- 
bers of any cross section whatsoever. These separate steps would be re- 
quired to analyze members of any particular cross section: 


1. Establishment of an expression or expressions relating the unit angle 
change to the applied moment. 

2. Integration of the unit angle change for various moment variations to ar- 
rive at expressions for total angle change and tangential deviation. 

3. Drawing up of curves relating deformations to applied loads. 

4. Once the curves of Step 3 have been drawn up, the remaining analysis 
would follow strictly the example problems presented. 


Several points should be mentioned in connection with this procedure. The 
rigorous treatment of I-beams would result in considerable difficulty. There 
are actually two distinct inelastic stages in an I-beam in bending: Under in- 
creasing moment, once the elastic capacity of the beam has been exceeded, 
first the elasto-plastic frontier will be in the flange, then in the web as the 
ultimate moment is approached. While the relation between unit angle change 
and moment for the latter stage is simple, the integration of the expression 
for the first inelastic stage is quite involved. 

Also, for rigorous results, each I-beam shape would require a separate 
set of calculations and curves, obviously an excessive requirement. 

On the other hand, it should be noted that, due to the smaller inelastic 
range of I-beams, approximate methods which simplify the moment-curvature 
relationships beyond the proportional limit yield much better results for I- 
beams than for beams of rectangular cross section. 


CONCLUSION 


The proposed method offers an “exact” solution of inelastic deflections 
simpler than others heretofore devised. It might prove useful in establishing 
and comparing the validity of any approximate methods of finding inelastic 
deflections which have been or might be suggested. 


NOMENCLATURE 


Fraction of length of beam from End A to point of maximum moment. 
Fraction of length of beam from End B to point of maximum moment. 
Distance from neutral axis of beam to extreme fiber. 


Dimensionless coefficient pertaining to tangential deviation: 
L2 
ML 


Dimensionless coefficient pertaining to angle change: 6 = Co —F- i 
Modulus of Elasticity. 
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End moment on beam, divided by Myp- 
End moment on beam, divided by Myp- 
Moment of Inertia of cross section of beam. 
Length of beam. 
Bending moment on beam. 
Bending moment causing yield point stress in extreme fiber of beam. 
Concentrated load on structure. 
Yield point stress. 
Width of cross section of rectangular beam. 
Uniform load on structure. 
Distance along beam. 
Distance from neutral axis of beam to elasto-plastic frontier. 
Deflection of Point A of member, normal to axis of member. 


Deflection of Point A of member, calculated by traversing deflected 
shape of structure from Point B. 


Deviation of Point A of member, with respect to a tangent from 
Point B of the member. 


Total angle change between Points A and B of the Member. 
Slope of member at Point A. 

Strain in extreme fiber of beam. 

Yield point strain. 


Curvature (unit angle change) of beam. 
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LATERAL BUCKLING OF I-BEAMS SUBJECTED TO 
UNEQUAL END MOMENTS 


J. W. Clark,l and J. R. Jombock,2 Junior Members, ASCE 
(Proc. Paper 1291) 


SYNOPSIS 


An experimental investigation into elastic and inelastic lateral buckling of 
aluminum alloy I-beams subjected to unequal end moments is reported. An 
approximate method is presented for determining the inelastic buckling 
strength of I-beams with this type of loading. Suggestions are made by which 


nonuniform bending could be handled more accurately in present design 
specifications. 


INTRODUCTION 


The strength of beams as determined by lateral buckling in the elastic 
stress range has been calculated for a wide variety of conditions of loading 
and support. Knowledge of the elastic buckling strength is not very useful 
for design purposes, however, unless means are available for extending the 
elastic buckling solution into the inelastic stress range. An analytical solu- 
tion for the inelastic buckling strength of beams is difficult because the re- 
sistance of the material to bending and twisting can vary both through the 
depth of the beam and along its length. The problem can be studied more 
easily by means of tests. 

Until recently, investigations of the buckling strength of beams, either 
elastic or inelastic, had been confined largely to cases in which the bending 
moment diagram was symmetrical about the center of the span. Recently, 
the mathematical solution has been obtained for the elastic buckling strength 


Note: Discussion open until December 1, 1957. Paper 1291 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Scciety of Civil 
Engineers, Vol. 83, No. EM 3, July, 1957. 
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of beams subjected to unequal end moments.(1) The experimental investiga- 
tion described in this paper was carried out to study both elastic and inelas- 
tic buckling of beams subjected to unequal end moments. Tests were made in 
which the values of the ratio, r, of the moment at one end of the beam to the 
moment at the other end were 1.0, 0.5, 0, -0.5 and -1.0. The beams were re- 
strained at the ends against rotation in the horizontal plane, rotation about 
their longitudinal axes, and displacement in the horizontal plane. 


Description of Specimens 


The specimens were lengths of 2014-T6 aluminum alloy extruded I-section 
with nominal dimensions of 2-3/8-in. depth, l1-in. flange width and .094-in. 
thickness for both flanges and web. Average measured cross-sectional di- 
mensions and section properties are shown in Fig. 1. 

The material used for these specimens is a typical structural aluminum 
alloy. Average mechanical properties of the beams in the longitudinal direc- 
tion are shown in Table I. The compressive yield strengths, tensile yield 
strengths and tensile strengths of all specimens were within 2 per cent of the 
average values listed. A compressive stress-strain curve for strains up to 
about 7 per cent is shown in Fig. 2. The modulus of elasticity for 2014-T6 


aluminum alloys is 10,600,000 psi and the modulus of rigidity three-eighths 
as great. 


Description of Tests 


The beams were tested in pairs in a 40,000-lb. capacity Amsler testing 
machine with the auxiliary fittings illustrated in Fig. 3. The pairs of beams 
when assembled in the end fittings provided a stronger unit for handling prior 
to testing and provided more resistance to rotation of the end fittings ina 
horizontal plane during testing than individual beams would have provided. 
The end fittings, shown in detail in Fig. 4, were designed to give an end con- 
dition approaching complete fixity. Counterbalances were used so that no 
moment was introduced into the specimens from the weight of the end fittings. 
The loading beams and knife-edge arrangement by which the end moments 
were applied to the specimens were designed to give the maximum practical 
resistance to horizontal rotation of the end fittings during the tests. 

Special care was used in assembling the two-specimen units and in the 
alignment of the units in the testing machine to prevent initial out-of- 
straightness of the specimens. Before testing the specimens, a load of 200 
lbs. was applied to each specimen and the out-of-straightness of the top 
flange (that is, the horizontal deviation from a straight line through the ends) 
was measured. For specimens failing by inelastic buckling, the maximum 
out-of-straightness for the 200-lb load was less than L/ 1800, L being the 
specimen length. For specimens failing by elastic lateral buckling, the maxi- 
mum out-of-straightness was found to be less than L/ 800 for the 200-lb load. 

About a third of the specimens tested had type A5-1 SR-4 electrical re- 
sistance wire strain gages mounted on the flanges at each end, top and bottom, 
at a distance of 1 in. from the end fitting. Clamped to these same specimens 
at the approximate location of maximum horizontal deflection were light 
aluminum alloy plates from which linen lines were led to 100-ohm Helipot 
continuous-rotation linear potentiometers. This arrangement made it possible 
to measure the lateral deflection and rotation of the specimens. The electri- 
cal signals from the strain gages, potentiometers, and a Baldwin SR-4 fluid 
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pressure cell in the testing machine hydraulic system were recorded with a 
Hathaway S8-C recording oscillograph, providing continuous records of de- 
formations of the specimens and the corresponding loads. 

For all specimens, the load was increased continuously until the ultimate 
load was reached. Fig. 5 is a photograph of one of the specimens taken after 
failure. 


Results and Discussion 


A. Elastic Lateral Buckling 


Table II lists the maximum applied bending moments and the correspond- 
ing nominal critical bending stresses (hereafter termed experimental critical 
stresses) at the centroids of the flanges, computed by the usual flexure formu- 
la. For moments greater than about 12,500 in.-lb., the experimental critical 
stresses are higher than the proportional limit of the material and, therefore, 
are slightly greater than the actual flange stresses at failure. The experi- 
mental critical stresses are shown graphically in Fig. 6, which also shows 
the corresponding theoretical critical stresses for beams with ends fixed 
against rotation in the horizontal plane.(1) The expression from Ref. 1 used 
for evaluating the critical moments from which the stresses were found is 


aL 


= maximum bending moment in beam when buckling occurs, 
in.-lb., 


Ur = a magnification factor reflecting the influence of the end 


moment ratio r and the end conditions, * i , 

Ko.1 = \ / 1+ atl = a factor dependent on the section proper- i 

(wL) ties, beam length, and end conditions, if 

2 

a2 = ree = an index of the section’s strength against warping, a 
in.2, 

distance between centroids of flanges, in., as 

= Ely = minimum flexural rigidity of the beam, lb. -in.2, 4 

E = modulus of elasticity, psi, q 


= minimum moment of inertia of the beam, in.4, 
= GJ = torsional rigidity of the beam, lb. -in.2, 
= modulus of rigidity, psi, 


= torsion constant, in.4, 


rao nas 


= laterally unsupported length of beam, in., 


One of the significant contributions of Professor Salvadori’s paper! is the 

demonstration that the effect of the conditions of restraint at the ends on the 

value of “, is sufficiently small that for design purposes My may be con- 
sidered to be independent of the end conditions. 
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a = reduced length factor (@ = 1.0 when ends are simply sup- 
ported in lateral direction and 0.5 when ends are fixed in 
lateral direction). 


In Fig. 6 the experimental critical stresses generally fall below the the- 
oretical curves. For the beams that failed at stresses below the proportional 
limit of the material, this difference presumably resulted from the fact that 
the ends of the beams were not completely restrained against rotation in the 
horizontal plane. This behavior has been found in previous investigations of 
lateral buckling.(2, 3, 4) An evaluation of this effect was made for all speci- 
mens of the present investigation that buckled in the elastic range. 

Fig. 6 was used to determine the experimental reduced length factor, @, 
for these specimens. For any specimen, a@ is the factor which, when applied 
to the actual specimen length L, gives a reduced length aL that may be used 
in Equation 1 to obtain the experimental critical stress. 

Fig. 7 is a plot of the experimental reduced length factors as a function of 
the actual beam length for all specimens that failed by elastic lateral buckling. 
A straight line was fitted to the plotted points by the method of least squares. 
The equation of this line was found to be 


« = .610 - .00134L (2) 
For beams of about 82-in. length, the above expression gives an experimental 
reduced length factor of 0.5 corresponding to complete fixity. Equation 2 does 
not apply to beams whose length exceeds 82 in. The fact that @ approaches 
0.5 as the length of the beam increases is to be expected since for very long 
beams the stiffness of the test apparatus relative to the stiffness of the beams 
is very great. For the shortest specimens buckling elastically, this expres- 
sion gives experimental reduced length factors ranging from about 0.55 to 


0.58, corresponding to lengths of 48 in. for r = -1.0 and 21 in. for r = 1.0, 
respectively. 


B. Inelastic Lateral Buckling 


The experimental critical stresses for almost half of the specimens were 
in the inelastic range. The failure in these cases occurred either by inelastic 
lateral buckling, by local buckling, or a simultaneous development of both. 
The type of failure was determined by visual observation. A comparison be- 
tween the test results and calculated critical stresses for inelastic buckling 
is illustrated in Fig. 8. 

To take account of the reduced stiffness of the beam in the inelastic stress 
range, Equation 1 can be modified by substituting a reduced modulus E' for 
the modulus of elasticity E. The result can be expressed 


For comparison with test results, it is convenient to rewrite Equation 3 as 
follows: 


= E= BC 
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where Gor stress at centroid of flange corresponding to Mey, psi, 


I moment of inertia of beam about principal axis normal to 
web, in. 


The quantity in brackets in Equation 4 may be thought of as the value of flange 
Strain at which elastic buckling would take place if the material behaved 
elastically up to failure. This quantity is plotted as the abscissa in Fig. 8 and 
critical stress as the ordinate. Values of @ were calculated from Equation 2, 
which was assumed to apply in the inelastic as well as the elastic stress 
range. 

Three different values of the reduced modulus E' were used in calculating 
the curves in Fig. 8. In each case, the value of the reduced modulus was as- 
sumed to be governed by the nominal stress at the mean fiber of the com- 
pression flange at the point of maximum bending moment in the beam. The 
lowest curve in Fig. 8 results from using the tangent modulus E; for E’. Use 
of the tangent modulus in the buckling formula gives good agreement with 
tests in the case of inelastic buckling of aluminum alloy columns,(5) and the 
tangent modulus curve has also been used in establishing design stresses for 
aluminum alloy beams. (6) 


The middle curve in Fig. 8 was based on the double modulus, Ep, given by 
the formula 


(5) 
(VE+ VB)? 


The double modulus curve was found to agree closely with test results ina 
previous investigation of lateral buckling of beams with uniform bending mo- 
ment.(3) 


The highest curve in Fig. 8 was derived by assuming E' to equal the 
secant modulus Eg. 

The curves in Fig. 8 are cut off at a stress of 70,100 psi, which is the cal- 
culated stress for local buckling of the compression flange for the test beams. 
In calculating the local buckling stress for the flange, the attachment of the 
flange to the web was considered to be about half way between the fixed and 
hinged conditions, and the effective modulus was assumed to be the secant 


modulus (Ref. 7, p. 353). The resulting approximate formula for the flange 
buckling stress is 


where tf thickness of flange, in., 
bf width of outstanding portion of flange, in. 


For these beams, Equation 6 gives a buckling stress equal to 0.024 Es, which 
is equivalent to 70,100 psi. One of the assumptions on which this value is 
based is that the stress is uniform along the flange. The flange buckling 
stress would be expected to be higher when the beam is subjected to nonuni- 
form bending moient, so that the maximum stress occurs at only one point. 
The nominal stress values at which failures occurred that involved local 
buckling, including cases in which lateral buckling and local buckling occurred 
simultaneously, varied from 69,200 to 81,600 psi. 
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Measured and calculated critical stresses are further compared by the bar 
graphs of Fig. 9, which show the ranges of the ratio of experimental to calcu- 
lated critical stress. The four bars for each value of end moment ratio in 
Fig. 9 correspond to the four sets of calculated critical stresses. 

The first set of calculated critical stresses was based on the assumption 
that the critical stress for any beam subjected to unequal end moments may 
be approximated by the calculated critical stress for the beam subjected to 
uniform bending, using the tangent modulus as the reduced modulus in calcu- 
lating inelastic critical stresses. The reason for comparing calculated criti- 
cal stress based on uniform bending with experimental critical stresses for 
beams with nonuniform bending moment is that design specifications(6, 9) for 
allowable stresses in beams are generally based on the uniform-bending case 
(or a loading condition that gives nearly the same buckling stress, such as 
uniform load on the top flange). Fig. 9 shows that the beams with r = 0 car- 
ried up to 96 per cent more bending moment (considering both elastic and 
inelastic buckling) than would the same beams subjected to uniform bending. 
When r was -0.5 and -1.0, bending moments at failure were as much as 160 
per cent greater than would be expected from considerations of the uniform 
bending case. This comparison shows that current design requirements for 
beams may be very conservative when the beams are subjected to reverse 
bending. 

The second set of calculated critical stresses was obtained from the elas- 
tic buckling solution for the type of bending moment that the beams were ac- 
tually subjected to, with the elastic modulus replaced by the tangent modulus. 
The comparison between the values obtained in this way and the experimental 
critical stresses is the same as the comparison between the tangent modulus 
curve and the test points in Fig. 8. It is reasonable to expect the tangent 
modulus curve to be conservative when applied to beams (Ref. 7, p. 165), and 
this is borne out by the test results, as shown in Fig. 9. The experimental 
critical stresses ranged from 8 per cent less to 39 per cent more than the 
calculated critical stresses. 

This second method of calculating critical stresses is essentially the same 
as the “equivalent slenderness ratio” method. The equivalent slenderness 
ratio for a beam is equal to the slenderness ratio of a column that would have 
the same elastic buckling stress as the beam. The equivalent slenderness 
ratio is used in conjunction with a curve of column buckling stress to deter- 
mine the buckling stress (either elastic or inelastic) for the beam. If the 
column curve is based on the tangent modulus, the resulting calculated beam 
strengths are the same as those determined by substituting the tangent modu- 
lus for the modulus of elasticity in the elastic buckling formula for beams. 
Similar results would also be achieved by applying the equivalent slenderness 
ratio method in combination with the straight-line formula for aluminum alloy 
columns,(10) since the straight-line formula agrees closely with the tangent- 
modulus curve. 

A third set of calculated stresses was obtained by the use of the double 
modulus in the elastic buckling solution. As in a previous investigation,(3) 
use of the double modulus gives good agreement between experimental and 
calculated critical stresses for beams under uniform bending moment, the 
experimental values being from 8 per cent less to 9 per cent more than the 
calculated critical stresses (see Fig. 9). For beams with nonuniform bending 
moment, the double modulus curve is conservative, as has been predicated, (8) 
since the modulus is reduced by the maximum amount in but a small portion 
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of the beam, while the remainder of the beam is stiffer. Including all the test 
results, the experimental critical stresses were from 8 per cent less to 34 
per cent more than the calculated values. 

Finally, a fourth set of calculated critical stresses was obtained by the use 
of the secant modulus in the elastic buckling solution. The computed values 
were cut off at the calculated local buckling stress. It can be seen in Fig. 9 
that the experimental critical stresses ranged from 12 per cent below to 19 
per cent above the computed values. Considering only the beams for which 
the end moment ratio was zero or less, the experimental buckling stresses 
are from 7 per cent less to 19 per cent above the values computed on the 
basis of the secant modulus. A rational basis for expecting the secant- 
modulus buckling stress to give a reasonable approximation to the true value 
for values of r between zero and -1.0 is explained in the appendix. 

Based on the results of these tests as summarized in Fig. 9, the critical 
stress for beams subjected to more or less uniform bending (r = 1.0 to 0.5) 
may be calculated by the use of the double modulus in Equation 4 and critical 
stress for beams with smaller end-moment ratios may be calculated by the 
use of the secant modulus in Equation 4. 


C. Recommendations Concerning Design Methods 


As the tests show, design specifications for beams based on the assump- 
tion of a uniform bending moment or a uniform top flange loading can be very 
conservative when applied to beams with unequal end moments. However, the 
extra strength of beams subjected to unequal end moments could be taken into 
account approximately without making any radical changes in current speci- 
fications. Present design specifications give formulas or charts that relate 
the allowable beam stress to a function of the dimensions of the beam, such 
as L/y B,/Sc inthe ASCE specifications for aluminum alloy structures(6) or 
ld/bt in the AISC specifications for steel structures.(9) These same formulas 
or charts could be applied to beams with unequal end moments if the function 
of the dimensions were multiplied by a coefficient whose value would depend 
on the ratio of the end moments, r. For example, the quantity L/y B1/Sc in 
the ASCE specifications would be replaced by CL/y B1/Sc. Values of the co- 
efficient C corresponding to various values of r would be given in a table or 
chart. 

The coefficient C described in the preceding paragraph is equal to 1A uy.* 
Fig. 10 shows minimum safe values of uy plotted versus r for beams with 
any condition of lateral support at the ends, as given by Salvadori.(1) It has 
been shown in the present investigation that these values of My apply approxi- 
mately to inelastic buckling, as well as to the case of elastic buckling for 
which they were calculated. The test results also help to show that an appre- 
ciable improvement in accuracy of design specifications would result from 
the introduction of a coefficient such as that suggested. 


CONC LUSIONS 


The following conclusions have been drawn from the results of this 


* If a coefficient of this kind were applied to the factor ld/bt used in the ASCE 
specifications, the value of the coefficient would be proportional to 1/u,. 
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investigation of the lateral buckling of aluminum alioy I-beams subjected to 
unequal end moments: 


1. The experimentally determined elastic buckling stresses are generally 
lower than theoretical critical stresses calculated on the basis of the assump- 
tion that the beam ends are fixed against rotation in the horizontal plane. 
Satisfactory agreement between experimental elastic buckling stresses and 
theoretical values is obtained by the use of an empirical reduced length factor 


in the theoretical solution. The reduced length factor varied from about 0.50 
to 0.58. 


2. Beams under unequal end moments can develop much higher stresses 
(as much as 160 per cent higher in these tests) than the same beams when 
subjected to uniform bending moment. Therefore, design specifications 
based on the case of uniform bending moment, or some other loading giving 
similar values of critical stress, can be very conservative when applied to 
beams with unequal end moments. 


3. Critical stresses in the inelastic range for beams under unequal end 
moments can be estimated conservatively by substituting the tangent modulus 
corresponding to the maximum stress for the modulus of elasticity in the 
elastic buckling formula. This method of calculation gives the same result as 
the use of the “equivalent slenderness ratio” method in conjunction with a 
curve of column buckling strength. Experimental critical stresses were 
from 8 per cent below to 39 per cent above the values calculated by this 
method. 


4. The use of the double modulus for the effective modulus in the elastic 
buckling solution gives inelastic critical stresses in reasonable agreement 
with test results for beams subjected to more or less uniform bending 

(r = 1.0 and r = 0.5) but results in conservative inelastic critical stresses for 
beams with smaller end-moment ratios. For all tests, the experimental 
critical stresses varied from 8 per cent less to 34 per cent more than the 
values calculated using the double modulus. 


| 


5. The use of the secant modulus for the effective modulus in the elastic 
buckling solution gives inelastic critical stresses in reasonable agreement 
with test results for beams with unequal end moments. Critical stresses for 
the test beams with r =0, -0.5, or -1.0 were from 7 per cent below to 19 per 
cent above the secant-modulus calculated values. 


= 


6. Current design specifications could be modified easily to take approxi- 
mate account of the greater strength of beams subjected to unequal end mo- 
ments as compared to beams loaded by equal end moments or uniformly dis- 
tributed load. The charts and formulas for beam strength already in use 
would apply to the case of unequal end moments if the function of the beam 
dimensions (for example, L/y Bi/Sc in Ref. 6 or ld/bt in Ref. 9) were multi- 
plied by coefficient whose value would depend on the end moment ratio, r. 
Values of the coefficient could be presented in a table or chart. 
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APPENDIX 


Application of the secant modulus to the calculation of buckling stresses 
for beams with nonuniform bending moment can be given some rational justi- 
fication by using the following reasoning: Let Ex be the effective modulus for 
the beam cross section at the point x. Assume that the buckling strength can 
be calculated by substituting an average value of 1/Ex for the quantity 1/E in 
the elastic buckling formula. To understand the meaning of this assumption, 
consider the expression for the strain energy of lateral bending, Vp, of the 
beam, including a variable modulus Ex. This expression is 


a2 7 


where ly is the lateral moment of inertia, u is the lateral deflection, and x is 
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distance along the longitudinal axis. An approximate value of Vp can be ob- 
tained by using the following expression instead of Equation 7: 


Rough calculations indicate that the use of Equation 8 rather than Equation 
7 can introduce fairly large errors in the values of Vp (as much as 20 or 30 
per cent), but the resulting errors in the value of critical stress are con- 
siderably smaller. Values of Vp given by Equation 8 can be either higher or 
lower than the correct values, depending on the variation of du/dx and Ex 
along the beam. 

If a modification similar to that described above is made in the expres- 
sions for strain energy of twisting and flange bending, the resulting equation 
for total energy of the beam is the same as the equation for elastic strain 
energy, (1) except that the quantity 1/E is replaced by the quantity in the de- 
nominator on the right-hand side of Equation 8. This quantity is the average 
value of 1/E; for the beam. Therefore, the resulting critical bending moment 
is equal to that obtained from the elastic solution by substituting the average 
value of 1/E, for 1/E. 

Assume that the effective modulus E, is the tangent modulus Et, corres- 
ponding to the stress in the compressive flange at the point x. Since the beam 
twists as well as bends when buckling occurs, the actual effective modulus 
should be greater than the tangent modulus. Thus, the use of the tangent 
modulus will tend to make the resulting calculated critical moments conserv- 
ative. The reciprocal of the tangent modulus is 


. 


do 


(9) 


where € = unit strain 


and C = unit stress. 


For a beam of length L, the average value of 1/Etx, designated by 1/Et » 


L (10) 
(te 


Assume that the flange stress varies linearly from 9, at one end of the 
beam to @] at the other end. For the beams of this investigation, the bending 
moment varied linearly. The assumption that the flange stress has a linear 


variation is only an approximation when some of the stresses exceed the pro- 
portional limit. Then 
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de. dedx_L_ de 


do dx do dx 


Combining Equation 10 and 11 gives 


“O 


Et L 


Integrating Equation 12 gives 


(13) 


Thus the average value of the reciprocal of the tangent modulus is equal 
to the reciprocal of the slope of the line connecting the two points on the 
stress-strain curve corresponding to the stresses at the two ends of the 
beam. If the end moment ratio is 0.0 or -1.0, the slope of this line is the 
secant modulus corresponding to the maximum stress in the beam. There- 
fore, when r has the value 0.0 or -1.0, and when the expressions for the vari- 
ous types of strain energy in the beam are assumed to be similar to Equation 
8, with E, being the tangent modulus at the point x, the resulting critical 
stress values are the same as though the secant modulus corresponding to the 
maximum stress were substituted for the modulus of elasticity in the elastic 
buckling formula. This does not mean that a theoretically correct solution of 
the buckling problem is obtained by substituting the secant modulus for the 
elastic modulus in the elastic buckling formula. It indicates only that there 


is reason to expect this use of the secant modulus to give an approximation 
of the correct solution. 


Ale 
icy 
4 
(11) 
de 
Me 
aa 
é,-€ 
a 
i 
\ 
a 
ii, q 
a 
: 
bas 
ag : 
q 
“4 
q 


1291-12 July, 1957 


TABLE I 


MECHANICAL PROPERTIES OF 2014-T6 
ALUMINUM ALLOY EXTRUDED I-BEAMS 


Tensile Strength,* psi 63,800 
Tensile Yield Strength,* 0.2 per cent offset, psi 60,300 
Compressive Yield Strength,** 0.2 per cent offset, psi 60,800 


Elongation,* per cent in 2-in. gage length 11.5 


Determined from standerd 1/2-in. wide rectengular tension 
specimens cut from web and machined to uniform thickness, see 
ASTM Standards, 1955, Part 2, HS8-54T, p. 1246-1260. 


Determined from 5/8-in. wide 2.63-in. long rectangular 
compressive specimens, cut from flenge and machined to 
uniform thickness (tested in subpress with Montgomery-Templin 
type jig), ibid., E9-52T, p. 1261-1267. 
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TABLE II 

| SUMMARY OF EXPERIMENTAL CRITICAL STRESSES 

| End Maximum Experimental 

Moment Applied Critical 

E Spec. Length, Ratio, Moment Stress, 
No. in. r in.-1b. | psi 

21 

0.5 60 1.0 2 800 10 600 
1 48 1.0 4 140 15 600 
1.5 42 1.0 4 700 17 600 
2 36 1.0 5 910 22 200 
2.5 30 1.0 7 530 28 300 
3 24 1.0 9 690 36 400 
3.5 21 1.0 11 820 44 400 
4 18 1.0 13 560 51 000 
5 12 1.0 15 270 57 
5.5 6 1.0 18 760 70 5008 
7 57 0.5 3 800 14 300 
” 39 0.5 6 740 25 300 
- 30 0.5 9 960 37 400 
9 21 0.5 15 080 56 700 
10 12 0.5 18 400 69 2008 
10.5 6 0.5 21 200 79 800 b 
ll 54 0 5 790 21 800 
12 45 0 8 120 30 500 
13 36 0 10 480 39 400 
14 30 0 13 050 49 100 
15 24 0 14 980 56 300 
16 18 0 18 430 69 300 
16.5 ) 0 19 670 74 000 
1” 72 -0.5 6 220 23 400 
18 54 -0.5 8 350 31 400 
19 36 -0.5 13 470 50 600 
20 27 -0.5 17 000 63 900 
21 18 -0.5 21 000 79 0008 
21.5 g -0.5 21 700 81 6006 
22 84 -1.0 4 710 17 700 
23 66 -1.0 6 520 24 500 
24 48 -1.0 9 310 35 000 
25 36 -1.0 13 000 48 900 
26 30 -1.0 14 940 56 200 
27 24 -1.0 17 540 65 900 
28 18 -1.0 20 100 75 700d 
29 12 -1.0 20 500 77 2006 


Failed by combined local buckling o ange and web and lateral buckling Bi 
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AVERAGE MEASURED 
OIMENSION DIMENSIONS, INCHES 


b 0.989 
d 2.361 
h 2.272 
te 0.0866 
ty 0.0888 


SECTION PROPERTIES BASED ON AVERAGE MEASURED DIMENSIONS: 
MOMENT OF INERTIA ABOUT X-X AXIS, Ty, 0.302 IN“ 
MOMENT OF INERTIA ABOUT Y-Y AXIS, ly, 0.0144 IN“ 
TORSION CONSTANT, J, 0.000955 in* 
TORSION BENDING CONSTANT, Cs, 0.0186 


DIMENSIONS AND SECTION PROPERTIES 
OF 2014-T6 ALUMINUM ALLOY EXTRUDED 


Ib. per sq. in. 


SPECIMENS CUT LONGITUDINALLY 
FROM FLANGES OF 2014-T6 
ALUMINUM ALLOY EXTRUDED I-BEAMS 


COMPRESSIVE STRESS, 


STRAIN, in. per in. 3 
AVERAGE COMPRESSIVE STRESS-STRAIN CURVE 
FOR ALUMINUM ALLOY EXTRUDED |-BEAM 
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FIG. 3 ARRANGEMENT FOR TESTING PAIR OF I-BEAM SPECIMENS 


(KNIFE -EOGES ARRANGED AS SHOWN ABOVE GIVE UNIFORM BENDING MOMENT, 
ret; FOR OTHER VALUES OF PLACE KNIFE-EDGES AS SHOWN BELOW) 
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FIG. 7 EXPERIMENTAL REDUCED LENGTH FACTOR 
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THE LATERAL RIGIDITY OF SUSPENSION BRIDGES 


I. K. Silverman,! A.M. ASCE 
(Proc. Paper 1292) 


SYNOPSIS 


The behaviour of a suspended structure under lateral loads is shown to de- 
pend upon the values of three dimensionless parameters 


For the important case of a uniform load over the entire span curves, appli- 
cable to moderately long span are presented for moment, shear, and deflec- 
tion coefficients as a function of the parameter A. The spring constant for 
small vibrations in the horizontal plane is also derived. 


Action under Static Lateral Forces 


When the suspended structure of a suspension bridge is deflected out of 
its plane by lateral forces, restoring forces due to existing suspender loads 
are brought into play. Since the suspender forces are a measure of the cable 
tension the behaviour of the suspended structure depends upon the cable stiff- 
ness in addition to its own stiffness. The lateral stiffness of a suspension 
bridge can be expressed in terms of a parameter which is the ratio of cable 
stiffness to lateral truss system stiffness. Due to this interaction, lateral 
loads acting on the suspended structure are transferred to the towers in part 
by the cable. 


The Uniform Distribution Assumption 


By making the simplifying, but inaccurate, assumption that the load trans- 
ferred to the cable is uniform, the condition for equality of deflection of the 


Note: Discussion open until December 1, 1957. Paper 1292 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 3, July, 1957. 
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cable and truss at the bridge mid span results in a direct solution of a load 
distribution between cable and truss systems. (See Appendix “A” for nota- 
tion.) 

Under the uniform distribution assumption and for uniform wind load over 
the full length of the ll which is the most important case, the restoring 
force & is given by [1] 


where A= 


Using and 102 


/ / 
z- 


Coefficients for deflection, bending moment and tower reaction at truss level 
are given by curves of Figures 5, 6, 7, and 8 marked “Uniform Distribution” 
for various values of A. 


The Elastic Distribution Method 


The Elastic Distribution Method as developed by Moissief and Lienhard [1] 
is a trial and error method for determining the distiibution of load between 
truss and cable. By assuming an initial distribution of the restoring forces 
z, the final distribution is obtained by preserving geometrical continuity be- 
tween truss and cable. This fundamental method has the disadvantage of all 
trial and error solutions in not being a direct solution. 


Differential Equation of Deflected Wind Truss 


Using the notation of Appendix “A” and Figures 1 and 2, the following re- 
lations may be written: 


fie” 

My= £A(t-X) 


2. Numbers in brackets refer to the Bibliography at the end of the paper—Ap- 
pendix B. 
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Eliminating 4 and u, the remaining unknown, v, is defined by the following 
differential equation: 


No solution of (1) in closed form is known. Power series, collocation, and 

finite difference methods [2], [3], [4] have been used successfully. None of 
these solutions however are in general terms which involve the parameters 
of the problem. 

An approximate solution of (1) in terms of a Fourier Series can be obtained 
using the Principle of Virtual Work, [5] [6]. For the uniform lateral load over 
the full span previously considered, an expression for v satisfying all bound- 
ary conditions can be written as 


LFA, SIN MAX AA, SIN 11. (2) 


Independent variations of the coefficients a, yields the following expression 
for a virtual displacement, du a 


du-= da, sia m1 + Sa, Si 17,1 


The Principle of Virtual Work requires that the following conditions be satis- 
fied: 


[L sia =0 

si (3) 


[Z SIT My LAL = O 


The expression for L under the integral sign is obtained by inserting (2) and 
its derivatives into (1). Unless (2) is a solution of (1) the resulting expres- 
sion will be different from zero. From (3) as many independent equations 
are obtained as unknown coefficients a,. Performing the integrations a set 
of simultaneous equations for ap are obtained, the general form of any equa- 
tion being given by: 


= if = (1) a 
ag 
ig 
1 
ty 


ASCE SILVERMAN 1292-5 


N= 


4 


i may have any odd integral value except n, depending upon the number of 


terms assumed in (2). For three terms in (2) the set of simultaneous equa- = 
tions obtained is the following: a 
4, 


(5a) 


byw, —-W/. Za (5c) 
ij 


The parameters in terms of which a, can be expressed by simultaneous 
solution of (4) are defined by 


Hhe . f= 


The right hand members of (4) are interms of @, B, W;, w c: For any wa Py 
given set of values of a, B, W;, We the simultencous solution of (5) gives ia 
numerical values of the coefficients apn. Thus when applied to Proposed es 
Design No. 2 as given in [1] with the following physical characteristics: 7? 

[= 4000 teet = 1150 /bs. per FF. 

h.= 509 feet wz = 150 /bs, per tt 

519 teet ur= 1300 pertt 
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the dimensionless parameters 


MO 0.1622 
22x ZOO0 


0.159! 


00 
L10,000x 4000. __ 

-/02 


a 


The numerical values of ay for a five term series (2) are: 


Kips/tt 
@,= 9.000 //7 /3/ 


Qs = 0.000 002, 24Z,0/ 
0.006 00/, O94, 3/ 
2.000,000, 398, / 


Moment and deflection curves are shown in Figures 3 and 4. 
The parameter of greatest influence is defined by 


A= 72 (6) 


For given values of a and B, for example o@ = 0,162 and B = 0.159, and 
with the relation Ww; = 0.9 w the simultaneous equations for a, reduce to 
those shown in Table I, where n = 1, 2, 5, 7 and 9 and A is taken as the 
variable. Solving these equations for various values of A the coefficients 

Ko, ki, K2, and ks are given in Figures 5, 6, 7, and 8. Table II gives a compil- 
ation of such coefficients for representative spans. Practically, the shape of 
the (A, k) curves are independent of a, B, W,, and Ww, values. For moder- 
ately long spans values of the coefficients from Figures 5 to 8, inclusive, are 
sufficiently accurate for design purposes regardless of the actual a, £, wt, 
and w. values, since actual values of o, 8, w;/W, are confined within 
narrow limits of variation. For A = 0 the (A, ” curves pass thru the 
points ko = 5/384; k, = 0.125; ke = 3/32; and ks = 0.5. For large values of 
X the values of k, may change sign from plus to minus. 


Oscillations in the Horizontal Plane 


Small vibrations in the horizontal plane about the position of equilibrium 
are described by means of the following differential equation: 
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(7) 


and m = mass of suspended structure per unit of length = p/g. Assuming 
modes of vibration of the form 


where W, = circular frequency 


the requirement that 
fF du =O 


leads to the following expression for 


+ 


2 
b 
- $ - Get 
C= Fatal weight ot cables per tt of bridge 


The spring constant defined as mw, becomes 
K= al xl 


The natural frequency is given by 


ht 
hie = 1.02 and 


n = 1 the natural frequencies for representative bridges are given in Table 
Ii. 


Using a value of X = 0.15 which corresponds to a value of 
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APPENDIX A 


coefficients of Fourier Series 

weight of all cables per foot of bridge 

acce.eration due to gravity 

sag of cable (Figure 1) 

sag of cable plus length of hanger at midspan (Figure 1) 
length of hanger at distance X from end of span (Figure 1) 


coefficient of deflection at span center (Figure 5) 


coefficient of moment at span center (Figure 6) 


Mo 


wy, 


coefficient of moment at quarter point (Figure 7) 
Mf 
wr, 


coefficient of tower reaction at truss level (Figure 8) 


— 


w, 


length of center suspended span (Figure 1) 


mass of suspended span per unit of length 


total vertical load carried by hangers per foot of bridge 

Wt - 

time variable 

deflection of cable at distance X from end of span (Figure 2) 


deflection of lateral wind truss at distance X from end of span 
(Figure 2) 


Wy + Wo 

horizontal wind pressure on cables per foot of bridge 

horizontal wind pressure on suspended structure per foot of bridge 
distance measured from end of span (Figure 1) 

vertical distance from top of tower to cable at distance (Figure 1) 

force of restitution due to inclination of hangers per foot of bridge 
modulus of elasticity = 29,000,000 lbs. per sq. in. 


| 
g = 
he 
ko = 
- Ely Ga 
= 
fo | 
| 
= 
= 
p = 
q 
t 
i 
u = | 
= 
w = 7 
i 


SILVERMAN 1292-17 


horizontal component of pull of all cables due to vertical loads. The 
increase in H due to wind loads is neglected. 


equivalent moment of inertia of lateral wind truss 
spring constant for horizontal oscillations 


= bending moment in wind truss, simply supported, due to wind loading 


dimensionless parameter = 
dimensionless parameter 
dimensionless parameter = 


dimensionless parameter = 


circular frequency 


he 3 n7r? 
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BUCKLING OF PLATES UNDER NON-HOMOGENEOUS STRESS 


P. P. Bijlaard,* M. ASCE 
(Proc. Paper 1293) 


SUMMARY 


Plastic buckling stresses are calculated for long plates, clamped at their 
unloaded edges and subject to longitudinal bending or eccentric compression 
in their plane. As in the author’s earlier work on simply supported plates to 
this end the governing partial differential equation is reduced to ordinary 
finite-difference equations, using second-order differences. The plastic 
buckling stresses of clamped and hinged flanges are calculated by the energy 
method. For clamped and hinged flanges also the elastic buckling stresses 
are calculated. Simple design formulas and graphs are presented for the 
plastic reduction factors with which the buckling stresses for the associated 
elastic cases have to be multiplied for obtaining the plastic buckling stresses. 


INTRODUCTION 


This paper is sequential to reference 1, where, after some discussion of 
the basic principles of the author’s theory of plastic buckling of plates,(2, 3, 4) 
the plastic buckling stresses were computed for long simply supported 

plates under eccentric compression or bending in their plane. The same 
method, using second-order finite difference equations, is used in the present 
paper for calculating the plastic buckling stresses of long clamped plates un- 
der the same loadings. It was shown in reference 1 that the resulting plastic 
reduction factor 7 is practically not influenced by the shape of the stress- 
strain curve and depends mainly on the maximum edge stress and edge strain. 
Therefore, as the most convenient way to calculate 7 , only material with a 
stress strain diagram as applies for mild steel was considered. (As to the 
actual applicability of the theory to mild steel and other metals the reader is 
referred to the “Conclusions.”) 

For computing the plastic buckling stresses of long clamped flanges an 


Note: Discussion open until December 1, 1957. Paper 1293 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 3, July, 1957. 
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enérgy method was used, expressing the variation of the buckling deflections 
in the lateral direction in terms of the normal modes of vibration of a 
clamped-free beam, as was done earlier for computing elastic buckling 
stresses in references 6 and 7. Since for the elastic range the present 
equations reduce to another form than used in reference 7 for calculating the 
elastic buckling stresses of clamped flanges, also the elastic buckling 
stresses were computed for this case. The plastic buckling stresses for 
hinged flanges were computed by the energy method as well as by a more 
Simple method, comparing deflecting and restraining moments with respect 
to the hinged edge, acting upon an elementary strip bdx, where b is the flange 
width and x is measured in the longitudinal direction. 


Long Clamped Plate Subjected to Longitudinal Bending or 
Eccentric Compression in Its Plane 


The Finite-Difference Buckling Equations 


Dividing the plate (Fig. 1) into imaginary longitudinal strips of width 


FIG. 1 CLAMPED PLATE SUBJECTED TO 
ECCENTRIC COMPRESSION IN ITS PLANE 


h = b/n, for each of these strips the plasticity parameters A, B, D and F can 
be considered to be constants, so that for each of the strips the differential 
equation(2, 3, 4) 


4 » 


applies. (1) E is the elastic modulus and I = t3/12, where t is the plate thick- 


ness. 0, is the appropriate compressive stress in the strip. Assuming the 
buckling deflection as 


w = Y sin ( ™/a) = Y sin Ax (2) 
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where Y is a function of y alone, and a is the half wave length, this differen- 
tial equation reduces to the following ordinary finite-difference equation for a 
pivotal point k (Equation (9) of reference 1): 


where 


12D, (Bek) 

-34D, - ah 

+ 30(Bt2F), ah + 

to, (5) 


Ox [Grp (6) 


Here 0xp is the maximum edge stress (Fig. 1). The plasticity parameters 
Ak, Bx, Dk and Fx in a pivotal point k are given by Eq. (22) and (23) of 
reference 2 or Eq. (21) and (21a) of reference 4. For pure compression, as 
occurs in all problems considered in the present paper 


B, [ E./E)]/¥ 
4/y 
| /( se) 
Ey/E) 
(E/E,)- 


where v is Poisson’s ratio. These values can also be read directly from 
charts in reference 5, as functions of the tangent and secant modulus ratios 
Et/E and Es/E, belonging to 0, at the pivotal point. However, as shown in 
reference 1, to take account of the variation in the Y-direction of P,, C4, 
Cox, and instead of their values ‘at the pivotal points, their statically 
equivalent concentrated values at these points were used, as determined by 
Eq. (12) and (13) of reference 1. 

With n = b/h = 5 (Fig. 1), assuming imaginary points outside the plate, 
Eq. (3) can be written down for the points 1, 2, 3, and 4 in Fig. 1, resulting 
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in the four linear homogeneous equations (14) of reference 1. For the present 
case with clamped edges the boundary conditions are (Fig. 1) 


Y 


5 6 


Inserting this into the four linear homogeneous equations they contain only 
the unknown deflections Yj, Yg, Y3 and Y4 and can yield deflections differing 
from zero only if the determinant of the system is zero. This leads to 


C,, 
“a 


ja 


Assuming a ratio a/b, the only unknown is this determinant is the product 
€ 252. With 


j 
TN (8) 
where N is the flexural rigidity of the plate, from Eq. (18) of reference 1, § 
[10 856 (9) 


Numerical Results 


To save space only the determinant for n = 5 was presented here. Actually, 
like in reference 1, € “b“ was solved from eighth order determinants, obtained 
by assuming n = b/h = 9. Expressing the strains €, (Fig. 1) as 


and assuming a stress-strain diagram like that of mild steel, ky was comput- 
ed for the case o@ = 2 (pure bending) with (E/E). = 0.6 at the upper edge (Fig. 
2b) and for o@ = 1, where 9, is zero at the lower edge, with (Es/E). = 0.5 at 
the upper edge (Fig. 2c). In both cases E;/E is zero at the upper edge. 
Choosing several half wave lengths, the minimum buckling stress coefficients 
ky and the optimum half wave length a were found. They are given in Table 1. 
the plastic reduction coefficients n are found as 7 = k,/ke, where ke is the 
buckling stress coefficient for the associated elastic case, that is the case 
with the same strain distribution from Eq. (10) as the plastic case. Hence for 
the associated elastic case 


6 6, [1-4/9] (11) 
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FIG. 2 CLAMPED PLATE SUBJECTED TO BENDING 
OR ECCENTRIC COMPRESSION IN ITS PLANE 


where (xe is the edge stress at the upper edge. The coefficients ke for a = 2 
and & = 1 from references 6 and 7 are 39.6 and 13.56, respectively. For the 
two cases in Table 1 this gives 7 = 22.3/39.6 = 0.563 and 7 = 5.95/13.56 = 
0.439, respectively, as also shown in Table 1. 


Table 1. Results of Calculations 


For @ = 0, that is, for pure compression, from Table 5 of reference 4, the 
plastic reduction factor is 


= 0.5987AD + 0-312(B+2F) 


Charts for the plastic reduction factors 7 9 for various cases of homogene- 
ous compression were made at the Bell Aircraft Corporation and were partly 
published in reference 9. Fig. 3 gives Ng from Eq. (12) for the present case 
of a clamped plate as presented in reference 5. Since for this case 1g is 
smaller than for a plate with simply supported edges, also, for example, for 
the case @ = 2 it was found that 7 is slightly smaller than (E,/E)e, while for 
the same case for a simply supported plate, in reference 1, 7 was found about 
equal to (Es/E)e. 

This applies for a diagram like that of mild steel. If, as another extreme, 
a hypothetical diagram like that of Fig. 4 is assumed, where E/E = Es/E, it 
may be seen from Eqs. (21) and (21a) of reference 4 or more directly from 
the charts of reference 5 that A, B + 2F and D become a little more than Es/E 
times their elastic values 1/(1- v2), Hence for any loading and boundary con- 
ditions the plastic reduction factor is a little more than Es/E. Therefore in 
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FIG. HYPOTHETICAL STRESS- 
STRAIN DIAGRAM WHERE 


E,/E = E,/E 


the present case, for all values of a, 1 would be somewhat more than E,/E. 


Also 11g crom Eq. (12) would be somewhat more than E;/E, as is seen direct- 
ly from Fig. 3. 


Design Formula 


To take all these circumstances into account the plastic reduction factor 
for any stress-strain diagram and values & between zero and two may be cal- 
culated sufficiently accurately from the formula 


= ke = 0.46a(E,/E) + (1-0-460)(%,), (13) 


(E,/E)e and (7 0)e refer to the upper edge of the plate, where the compressive 
stress is highest. "9 is given by Eq. (12) and can be read from Fig. 3 for 
Es/E and E;/E at the upper edge. 7 from Eq. (13) is given in Table 1 as n 
(13), for comparison with the computed 7 values. For @ = 0, for any stress 
strain curve, Eq. (13) yields 7 = (NQ)e, as given by Eq. (12). 

The plastic buckling stress coefficient kp is equal to 7 ke, where Ke refers 
to the associated elastic case, having the same strain distribution (@) as the 
actual plastic case. For convenience Ke from references 10 and 11 is given 
graphically in Fig. 5. 


Long Clamped Flange Subjected to Longitudinal Bending or Eccentric 
Compression in Its Plane 


Derivation of Buckling Condition 


This case is most conveniently solved by an energy method, expressing the 
variation of the deflection in the Y-direction (Fig. 6) in the normal modes of 
vibration of a clamped-free beam, as was done in reference 6 for clamped 
plates and in reference 7 for clamped flanges, both for elastic buckling. 
Hence the buckling deflection 
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5 ELASTIC BUCKLING STRESS COEFFICIENT k, FOR A CLAMPED 
PLATE UNDER VARIOUS LOADINGS 


STRESSES BEFORE 
BUCKLING AND DEFLECTION 
AT INCIPIENT BUCKLING 
FOR CLAMPED 

FLANGE 
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X = sin(t/a) = sin AX (15) 


where a is the half wave length. From references 7 and 12, imposing the con- 
dition that 


2 
where b is the flange width (Fig. 6), 


> cos y/b) - cosh(p. - [since y/s) - sinh(p, (17) 
with 
Cosp, + coshp 
sin B, t+ sinh 
Here Pp has to satisfy the condition 


As was expected, also in the plastic range this series converges very satis- 
factory, so that only three terms are needed in Eq. (14). 


In the elastic range, from the principle of virtual displacements, at in- 
cipient buckling the energy 


U < constant + V, - A (20) 


reaches an extreme value. If in here Vj and Ve are the work done during 
buckling by the internal stresses and by the external forces, respectively, 
this remains true for plastic buckling. This is evident if one considers that 
the relative decrease of the buckling stress in the plastic range is due to the 
decrease of the secant and tangent moduli, but is not influenced by the fact 
that the plastic deformations are not reversible, since nowhere reversal of 
stress occurs. In other words, the plate behaves as an elastic plate with the 
same decreased rigidities as the plastic one, so that the principle of virtual 
displacements remains valid. 


From Eq. (59) of reference 13 the work done by the internal stresses dur- 
ing buckling is 
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where A, B, D and F are the plasticity parameters mentioned in the foregoing. 
The external work during buckling is, considering compressive stresses as 
positive, 


V, = (ae) dxay (22) 


where o, is a function of y (Fig. 6). From Eq. (14), (20), (21) and (22), fora 


. plate length a 
+ DX (2 anYu) + #FX (2 am Yo 
2 
- (Zam dx dy (28) 
i where a prime and a dot mean differentiation with respect to x or y, respec- 
7 tively. Since U should reach an extreme value, aj, ag, etc. should be chosen 
such that 
/ 0a, = O (24) 
where i = 1, 2, etc. This leads to equations of the general form ; 
[flex + BXX( YE anh) i. 
+ DXY, Yn YE am Yo] 
ii 
t 6, X am Yon dxdy = 0 (25) 
Integration of this equation in the X-direction gives, after division by a/2 and i 
- denoting all term numbers of the sums by m again me 
ij 
ii 
+ DY Z 4 Fat Zam Yu ] dy 
a 
2 
-tr A dy = 0 (26) 
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Case 1 


First the case of pure bending of a flange in its plane is considered (Case 
1), assuming a stress-strain diagram as that of mild steel and a stress dis- 
tribution as shown in Fig. 7a, where at both edges Es = 0.5 E. Here in the 
plastic parts 0, =t Oxp and in the elastic parts 0, = Oxy) [-2+ (4y/b)], so 
that in Eq. (26) 


2 


4b 


The normal modes Y,, are orthogonal functions, so that, if i 4 m, 


(28) 


ay 


Dividing in Eq. (26) the expression to the left by A 2772 EI/b2 one obtains, us- 
ing Eqs. (8), (16), (27) and (28), 


dy - kp, [01-»") 0 

where V is Poisson’s ratio and 


4, b 
‘ ab 
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The first term of Eq. (29) can be more conviently written as 
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yes 


where the subscripts e and p refer to elastic and plastic, respectively. For 
elastic deformation, where Eg = Et = E, with v = 0.3: A = D = 1/(1 - v2) = 
1.099, B = v/(1 - v2) = 0.329 eye 1/(2 + 2v) = 0.385. 

Hence 
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% dy (33) 


t 0-77) am 1.0998 a S:] b 
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Jmi 2nd Rypj have been calculated in reference 7, from which, for m # i, 
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and for m =i 


hi - pC, ) (38) 


_ (39) 
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The expression Sj can be calculated by integration by parts. One obtains, 
since at y = 0, Y = Y'= 0 and at y = b, y" = Y""' = 0, 


b :: 
ay ay (m/e) % 4 


di 

+ 

(34) 
where 

6 

[= 
il 

i 
j 
ail 


BLJLAARD 1293-13 
Hence, in connection with Eq. (16) and (28), from Eq. (35), 


Using three terms of the series in Eq. (14), from Eqs. (18) and (19) 


Py = 1.85710 ; pg = 4.69474 ; pg = 7.85476 } 
= 0.73410 ; Cg = 1.01845 ; Cg = 0.999224 


and from Eq. (36) = (39) 


Ji 0.17391 ; J12 -2.38018 ; J = 5.56315; 

J94 0.37969 ; Joo -2.69482 ; Jo3 = -1.83184; 

J31 0.31703 ; Jgo 0.65440 ; J33 = -9.27126 (42) 

0.94184 ; = R21 = -1.49543 ; Ry3 = = 0.79871; 

R22 6.57043 ; = R32 = -4.53002 ; Rg3 = 15.63312 
Choosing values for §, the right hand side of Eq. (34) can then be calculated 
in terms of aj, ag and a3, for m = 1, 2 and 3 andi = 1, 2 and 3. The second 
and third integral to the right in Eq. (33), containing the plasticity parameters 


A, B, D and F, that vary with y, were evaluated by Simpson’s rule, dividing 
the plate width b into 8 equal spacings (Fig. 7a). A, B, D and F in the point 


Case |,a(z-8-y)=2 Case 2, %(z=6-y)=1 Gase 3, (z=y)=!1 
FIG. 7 CLAMPED FLANGE SUBJECTED TO BENDING OR ECCENTRIC COMPRESSION IN ITS PLANE 


0, 1, 2, 6, 7 and 8 follow from Eqs. (21) and (21a) of reference 4 or can be 
read from the charts in reference 5 as functions of E,/E and Et/E. In all 
these points E; = 0 and Es/E is 0.5 for points 0 and 8, 2/3 for points 1 and 7 
and 1 for points 2 and 6. The ordinates Yj and Ym in the points 0 to 8 follow 
from Eq. (17); their derivatives are obtained by simple differentiations. 

The integrals in Eq. (31) lead to rather extensive expressions, so that it 
was more convenient and accurate to evaluate them too by Simpson’s rule, 
especially since the Y,, values for points 3, 4 and 5 had to be calculated any- 
way for evaluating the first term in Eq. (29) for the next Cases 2 and 3. 

Hence Eq. (29), for chosen values of 8, leads to 3 homogeneous linear 
equations in the 3 unknown coefficients aj, ag and a3. Equating the denomina- 
tor determinant of these equations to zero, Kp is the only unknown and can be 
solved from a cubic equation (Approximation Il). The convergence of the 
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solution can be checked by assuming as a first and second approximation ag = 
ag = 0 and ag = 0, respectively (Approximations I and Il). 


Table 2. Values y°) for Case 1 


Table 2 shows that the convergence is very satisfactory. Kp = 0.91 (1.20) = 
1.092 at the optimum half wave length a = 1.45b. 
Case 2 

Similar calculation were made for eccentric compression with a stress 
distribution as shown in Fig. 7b, where at the free edge Es/E = 0.5 and Ey = 0. 
Eq. (26) also applies here, but now in the plastic range 0, = %yp and in the 
elastic range 0x = (2y/b) ?xp, So that in the last term of Eq. (28) 


b 
° ie 


Hence, in Eq. (29), where yj; remains given by Eq. (30), except that A, B, D 
and F have values that differ from those for Case 1. 


2 am Yoo] 4 (44) 


Further, instead of Eq. (33), 
(45) 
° a ane 
The first integral to the right is given by Eq. (34). All other integrals were ' 
calculated by Simpson’s rule. Eq. (29) yields the results given in Table 3. i 
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Table 3. Values k/(1-¥*) for Case 2 


Hence Kp = 0.91 (0.81) = 0.737 at the optimum half wave length of about a = 
1.11b. 


Case 3 


For this case the assumed stress distribution is shown in Fig. 7c. Here 
at the clamped edge Es/E = 0.5. In a similar way as in the other cases one 
finds that in Eq. (29), where yj remains given by Eq. (30), 


(46) 


Further, instead of Eq. (33), 


The results for this case are given in Table 4. 


Table 4. Values k,/(1- y*) for Case 3 


- 
Re 
aq 1.00 0.838 0.838 0.835 
1.10 0.828 0.828 0.812 
1.20 0.831 0.831 0.828 
4 / i ey 
y2) 
im 3.350 2.702 2.701 
1.7 3-305 2.695 2.695 
1.8 3.286 2.706 2.706 
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This gives Kp = 0.91(2.695) = 2.452 at the optimum half wave length a = 1.7b. 
Elastic Buckling Stresses for the Associated Cases 1a, 2a,and 3a 


In order to express the strain distribution in Cases 1, 2 and 3 by a uniform 
formula it is convenient to write 


= [1-(«2/b) (48) 


where €,, is the maximum edge strain, occurring at one of the edges, and z 
is the distance from that edge. Hence, from Fig. 7, for Cases 1 and 2, z = b-y 
and for Case 3, z= y. Further for Cases 1, 2 and 3, a@ is 2, 1 and 1, respec- 
tively. For the associated elastic Cases la, 2a and 3a 


\ 


6, = 6,[1-(«2/6)] (49) 


Using Eqs. (14) - (20) the elastic buckling stresses for Cases la, 2a and 3a 
were calculated in reference 7. However, these calculations are partly based 
on Eq. (19) of that reference, which, to the present writer’s opinion, does not 
apply for the present boundary conditions. Therefore these elastic buckling 
stresses were calculated by reducing the formulas, used in the foregoing for 
calculating the plastic buckling stresses, to the elastic range. For Case la, 
with z = b-y and @ = 2, 


6, = 6. “1 +(23/b)/ (50) 


and in Eq. (26), in connection with Eqs. (16) and (28), 
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From references 7, for m # i, 4 
and for m = i, 
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from which 


Hy, = 0.80654 ; Hy2 = Ha) = -0.15343 ; Hy3 = Hg, = 0.02032; 
Ho2 = 0.59412 ; H23 = Hgo = -0.19090 ; = 0.53239 


Eq. (29) now becomes 


where kK, is determined by 


k Bt (97) 


The first term in Eq. (56) is given by Eq. (34) and from Eqs. (26), (51) and 
(56) 


=(-a, +22 mi) b (58) 


Hence Eq. (56) can be evaluated. The results are given in Table 5. 


Table 5. Values k,/(1- v") for Case la 


This gives kg = 0.91 (2.362) = 2.149 at an optimum half wave length a = 1.66b. 
As may be easily verified, for Case 2a, with z = b-y and @= 1, in Eq. (56) 


bY am Hmi (59) 


which leads to the results given in Table 6. 


ASCE BILJLAARD 1293-17 
— 
% 0 (56) 
= 

I III 
1.57 2.378 26377 2.374 
1.6 2.372 2.371 2.369 
1.63 2.367 2.366 2.3640 
1.7 2.367 2.366 2.3641 


1293-18 July, 1957 


Table 6. Value k,/(1- y*) for Case 2a 


Hence Ke = 0.91 (1.778) = 1.618 for a = 1.66b. 
For Case 3a, with z = y and @ = 1, in Eq. (56) 


Q ay b (60) 


The results for Case 3a are given in Table 7. 


Table 7. Value k,/(1- y*) for Case 3a 


air 
4 
This gives kg = 0.91 (6.517) = 5.930 at a = 1.57b. Aer 
Design Formula 
The values kp and ke, calculated for Cases 1, 2 and 3 and their associated a 
elastic cases la, 2a and 3a, respectively, are assembled in Table 8. For all a 
cases, at the most loaded edge (E,/E), = 0.5 and E/E = 0 (Fig. 7). The re- 4 
sulting plastic reduction factors 7 = kp/Ke are given in the sixth column. * 
Also 7 for homogeneous compression (a@ = 0) is listed. It is calculated as the . 
ratio 0p/o¢ from Table 5 of reference 4, using Eq. (21) and (21a) of that a 
reference or can be read directly from a chart in reference 5, for Eg/E = 0.5 4 
and E,/E = 0. This chart is given here as Fig. 8. a 


1.6 1.803 1.784 1.7816 “a 
1.63 1.800 1.781 1.7789 
| 1.799 1.782 1.7798 
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FIG. 8 PLASTIC REDUCTION FACTORY, FOR A CLAMPED FLANGE 
UNDER HOMOGENEOUS COMPRESSION, AS A FUNCTION OF 
and 
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Table 8. Results of Computations 


For any value of o@ in Eq. (48) the reduction factor 7 with respect to the 
associated elastic case, having the same strain distribution as the plastic 
case, may be assumed to be given by the third degree curve passing through 
the four known values of 7 from Table 8. This gives 


= (%,). + (0-370 0-255 (61) 


where (Es/E). is the secant modulus ratio at the edge-where 0, has its high- 
est value and (75)e applies for homogeneous compression with a stress and 
strain equal to the actual maximum edge stress and strain, that is, belonging 
to Es/E and E;/E at the maximum loaded (in compression) edge. 1g can be 
read directly from Fig. 8. For a, the parameter of the strain distribution in 
Eq. (48), in Eq. (61) the plus sign has to be used if z = b-y, that is, if the 
maximum strain occurs at the free edge. If z = y, that is, if the maximum 
strain occurs at the clamped edge, a has to be inserted into Eq. (61) asa 
negative value. Actually Eq. (61) has been chosen such that, for ao = 2, 7 = 
(E,/E)e, instead of 2% more, as would follow from Table 8. 

Further, in Table 8, ke(7) and ke(11) are the elastic buckling coefficients 
Ke from references 7 and 11, respectively. The values in reference 11 were 
calculated from finite-difference equations based on the equation of the 
funicular polygon, as developed in references 14 and 10. Both values ke(7) 


and ke(11) were calculated with a Poisson’s ratio v= 0.3, as was also used in 


the present paper. From Table 8 the present values k, are in better agree- 
ment with k,(11) than with ke(7). 

For convenience Ke is given graphically in Fig. 9. No data is available be- 
tween ot(z = y) equal to zero and one, but in drawing the pertinent curve be- 


tween the known points it was taken into account that its continuation should go 
through k, = 0 for a(z = b-y) = 1, because there ke for the compressive stress 


at the clamped edge, where a(z = y) refers to, is zero. 


Hence kp = ke, where ke is given by Fig. 9 and 7 by Eq. (61). In the lat- 
ter 79 can be read from Fig. 8. 


Long Hinged Flange Subjected to Longitudinal Bending or 
Eccentric Compression in Its Plane 
Derivation of Explicit Formula for Plastic Reduction Factor 


Similar to the case of homogeneous compression (@ = 0) for a long flange 
it is sufficiently accurate to assume a buckling deflection 
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FIG. 10 HINGED FLANGE SUBJECTED TO BENDING OR ECCENTRIC COMPRESSION IN ITS PLANE 


(for Z=y) (for z=&-y) 
0.5 0.5 z 


FIG. 9 ELASTIC BUCKLING STRESS COEFFICIENT k, 
FOR A CLAMPED FLANGE UNDER VARIOUS 
LOADINGS 
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W Ky sin ( 


Using the condition that at incipient buckling the work Vj done by the internal 
stresses must be equal to the work Vg, done by the external forces, or 


(63) 


from Eqs. (21) and (62), for one half wave length 


| F ay (64) 


From Eq. (21) of reference 4 


F / 


Case 1 


Assuming pure bending and a stress strain diagram like that of soft steel 


(Case 1), with plastic zones of arbitrary width d, from reference 4 and Fig. 
10a, 


e = (E/E,) -! = Ee (66) 


where 


With v = 0.3, from Eq. (65) and (66) 


a/b 


(68) 
[ Fay, = b 


In the elastic part, where e = 0, F = 0.385, so that 


b 


n 43+ d 
+ 4 | b (69) 


From Eqs. (22) and (62) 
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From Fig. 10a, in here 


a L 


Inserting Vj and Ve into Eq. (63), it yields, with EI = (1 - v 2)N, 


6 i-ad) 3a5+(4 4)tn sted (72) 
bt 


130-26) 


For the associated elastic case, for which in Eq. (49) z = b-y and @ = 2 (pure 
bending), from reference 7, 


xe 


73 
(73) 


so that for this value of a 


Ore 13(7-2¢) 


According to Eq. (66) in here 6 can be expressed in the secant modulus ratio 
(Es/E), at the compressed edge, where y, =a: 


(6/4). | (75) 


The plastic reduction factor ng from Eq. (74) is given in Table 9 for various 
ratios (Es/E)e in the column marked z = b-y, @ = 2. 


Table 2: Plastic Reduction Factors 


In contrast with the preceding cases, for pure bending (@ = 2)” is here 
greater than (Es/E)e, especially for small ratios (Es/E),. In reference 1 it 


4 
4 
(71) 
% 
ee. 
4 
Gal 
q 
bat 
0.875 0.884 0.890 0.860 0.865 
0-750 0.768 0.785 0.722 0.705 14 
0.625 0.696 0.726 0.590 0.580 
0.500 0.605 0.645 0.465 
0375 0.505 0.550 0.342 0.320 
0.250 0.390 0.432 0.224 0.210 
= 
4 
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was shown, uSing its Fig. 13, that for a simply supported plate 7” will be about 
equal to (E,/E)e. As explained in the foregoing, for a clamped plate 7 is a 
little smaller. For a clamped flange, from Table 8, where for all cases 
(Es/E)g = 0.5, 7 is indeed about equal to (Es/E)e. 

For the present case of a hinged flange, however, the same reasoning is 
not valid. It is based on the assumption that the shape of the deflection for 
the plastic case is the same as for the elastic case, so that the restraining 
forces are reduced in the ratio of about (Es/E),. In the present case for both 
the plastic and the elastic case the deflection may be assumed as Eq. (62), but 
the deflecting forces depend on the fourth order partial derivatives of the de- 
flection, which, for an infinitely long plate, from Eq. (62) are all zero or in- 
finitely small of higher order. Hence, with the same approximate deflection, 
the restraining forces in the plastic case may still have a distribution differ- 
ent from (Es/E)e times those in the elastic case. By considering a strip bdx 
of the flange (Fig. 11) it is seen that its equilibrium requires that the moment 
Map, given by the deflecting forces D about the hinged edge y = 0, must be 
equal to the difference Mrp between the twisting moments at both edges x = 


const. of the strip. Using Eq. (62), the deflecting forces acting on an element 
dx dy of the strip are 


( = (w/a )t 6, Kysin(Tx/a)dxdy 
(76) 
C, y 4y 


where Cj, is a constant for the strip considered. Hence the total moment 
about the hinge is (Fig. 11) 


4 
Map at y 4y 


where the integral is given by Eq. (71). 
From Eq. (22) of reference 4 the difference between the twisting moments is 
(Fig. 11) 


b 
Mp =- (M,, dxdy = -2 =, [Fas (78) 


where the integral is given by Eq. (69). At incipient buckling 


M, =M 
r 


(79) 


Denoting the deflecting and restraining moments in the associated elastic case 
aS Mge and Myre, since also Mge = Mre, 


; 
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In a similar way as was calculated in Eqs. (71) and (77), in the associated 
elastic case (pure bending) 


b ‘ 
‘ 3 Z (81) 


where 0,, is the compressive edge stress. In the elastic range, with e = 0, 
from Eq. (65) F = 0.385, so that, using Eq. (78), 


M_./Co = 0.385 b (82) 


Hence in Eq. (80) Map is given by Eqs. (71) and (77), Mrp by Eqs. (69) and 
(78), Mge by Eq. (81) and Mre by Eq. (82). Insertion into Eq. (79) leads to the 
same expression for 0xp/dxe =" as given by Eq. (74). 

In the same way expressions for 7 were found for Cases 2 and 3. 


Case 2 
For Case 2, where in Eq. (48) z = b-y and @ = 1 (Fig. 10b) 


where now 


§ = d/b = 1-(Es/E)e (84) 
Case 3 
Here in Eq. (48) z = y and @ = 1 (Fig. 10c) and 


1301-0) 


where 


5 = d/b = 1-(E,/E)e (86) 


Values for ny from Eq. (83) and (84) are given in Table 9 under z = b-y, @ = 1. 
Those from Eq. (85) and (86) are shown under z = y, @ = 1. For homogeneous 
compression (a@ = 0), from Table 5 of reference 4, 79 = 2.6 F, where F is 
given by the present Eq. (65). In Table 9 values for ng are given under @ = 0. 
Using the values in Table 9 in Fig. 12 the reduction factor nq is plotted ver- 
sus a(z = b-y) and a(z = y) for various ratios (Es/E),. at the edge with highest 
compressive strain. Actually these ratios 1g are valid for the special stress- 
strain diagram assumed. 

As another extreme, let the stress-strain diagram be given by the straight 
line relation in Fig. 4, in which hypothetical case at any stress the strain con- 
sists of an elastic part € , and a plastic part € p=P €g; where p is a constant. 
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For any stress distribution (Fig. 13) then the ratio e = € /E, = p will be ¥ 
equal to that at the maximum compressed edge, so that from Eq. (65) every- = 
where F will be equal to Fe at that edge. Since in the elastic range F = 0.385, 
this makes the restraining moment Mpp in Eq. (78) Fe/0.385 = 2.6 Fp times 
that for an elastic flange with the same buckling deflection, so that in Eq. (80) 
Mrp/Mre = 2.6 Fe. 

The stress distribution for the plastic case is identical to that for the 
associated elastic case, so that the deflecting moment Mgp, as meant in Eq. 
(77), is Oxp/Oxe times that in the elastic case, 0xp and Oxe being the edge 
stresses in the plastic and elastic case, respectively. Hence from Eq. (80), 
for any value @ in Eq. (48), 


% = Sep ce = 426F, (87) 


where Fe refers to the most compressed edge. 


Design Formula 


Hence for any actual stress strain diagram, given by the solid curve OPQ 
in Fig. 14, it seems appropriate to obtain the plastic reduction factor for the 
maximum compressive edge stress 0xp as 


Here A, is the area between the actual stress-strain curve OPQ up to the 
maximum edge stress Oxp and the hypothetical diagram, given by the dashed 
straight line OQ and Ag is the area between the actual curve OPQ and the 
broken dashed line ORQ, which has the same shape as the diagram of soft 
steel. Further fy is given by Fig. 12 and F¢ is given by Eq. (65), where e is 
related to Eg/E by Eq. (66). Both Ng and Fe have to be determined for 
(Es/E)e at the most compressed edge. For a diagram like mild steel, Ag = 0, 
so that Eq. (88) yields 7 = Ng, as it should. For pure compression (a@ = 0), 
from Table 5 of reference 4, Ng =o = 2.6 F, so that from Eq. (88) for any 
stress strain diagram 7 = 2.6 F, in accordance with Table 5 of reference 4. 

The plastic buckling stress coefficient Kp is 1 Ke, where the buckling stress 
coefficient ke for the associated elastic case is given graphically in Fig. 15. 
As may be easily calculated in the same way as xp was calculated in the 
foregoing, for cases where the maximum compressive stress occurs at the 
free edge, that is, where in Eq. (49) z = b-y, 


ke (89) 
€ 2 
If the maximum compressive stress occurs at the hinged edge, so that in Eq. 


(49) Z=Yy, 


ke 46-8 (90) 
T*( 4-30) 
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xe ° & 


FIG. 13. HINGED FLANGE FIG. 14 ARBITRARY 
UNDER ECCENTRIC STRESS- 
COMPRESSION IN 
ITS PLANE, IF 
E,/E = £,/E (FIG. 4) 


(for z= t-y) 


FIG. 15 ELASTIC BUCKLING STRESS COEFFICIENT k, FOR A HINGED 
FLANGE UNDER VARIOUS LOADINGS 
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For ao = 1 and 2 Eqs. (89) and (90) yield values ke that practically agree with 
those calculated in reference 7. 


CONC LUSIONS 


As follows from experiments referred to in reference 1, for plates under 
homogeneous compression the author’s theory is in excellent agreement with 
tests on aluminum alloys. Also for other metals where Et is always positive, 
so that no yielding in layers occurs, the same good agreement may be ex- 
pected. For mild steel that agreement is satisfactory if the stresses are not 
at the lower yield point, that is, not in the horizontal part of the stress- 
strain diagram between upper yield point and the strain hardening range. In 
that region, where 9 = Pys, in the case of homogeneous compression, one can 
only speak of a critical strain, and from experiments reported in reference 15 
the actual critical strain is higher than predicted by the theory, that is, the 
theory is conservative. This is quite understandable, since it assumes con- 
tinuous plastic deformation, while at the lower yield point yielding in discrete 
layers occurs, leaving elastic blocks between them. This, of course, makes 
the material more rigid than would follow from the theory. This same phe- 
nomenon hampers the buckling of columns in this region, as was pointed out 
in reference 16. Hence the present results may be expected to be accurate 
for alloys of aluminum, magnesium, titanium, etc. They will also be accurate 
for steels that do not have a horizontal part in their stress strain diagram 
and for other steels if the stresses are nowhere at the lower yield point. If 
the stresses are all or partly at the lower yield point, the present results are 
conservative. 

For all loadings and boundary conditions the plastic buckling stress coeffi- 
cient in Eq. (8) is given as 


kp YM ke (91) 


For clamped plates and an assumed stress distribution, with maximum edge 
stress 0 xp , ke refers to the associated elastic case, for which the stress 
distribution is given by Eq. (11), where @ is determined by the strain distribu- 
tion in the plastic case, as given by Eq. (10). In both Eqs. (10) and (11) y is 
measured from the edge with the highest compressive stress. k, can be read 
from Fig. 5. The plastic reduction factor n is computed from Eq. (13), where 
(Es/E)e and (Ng)e refer to the edge where the maximum compressive stress 
is acting. 70 is read from Fig. 3 for Eg/E and E;/E at that edge. From Eq. 
(23) of reference 1 the required plate thickness t (without safety factor) is 


For a clamped flange or a hinged flange the stress distribution in the 
associated elastic case is given by Eq. (49), with the same @ value as occurs 
in Eq. (48) for the strain distribution in the plastic case, as explained in the 
text following that equation. For a clamped flange kp can be read from Fig. 9. 
The reduction factor 7 is given by Eq. (61), where (Es/E)e and (N0)e refer to 
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the edge with maximum compression. 79 can be read from Fig. 8 for Eg/E 
and E,/E at that edge. 

ke for a hinged flange follows from Fig. 15 and 7 is determined by Eq. (88) 
as explained in the text following that equation. 

For both clamped and hinged flanges the required thickness t is given by 
Eq. (92). 

In reference 17 the case of simply supported plates, under combined shear 
and eccentric compression in this plane, is dealt with. 


This investigation was carried out at Cornell University under the sponsor- 
ship and with the financial assistance of the Office of Ordnance Research of 
the U. S. Army. Most of the numerical work was done on a desk calculating 
machine by Mr. Mohammad Taheri or on I.B.M. machine by Miss. R. B. 
Lahayne, part of it both ways. 

The author wishes to express his thanks and appreciation for the valuable 
co-operation of all concerned. He also is indebted to the Bell Aircraft Cor- 
poration for its permission to present the graphs in Fig. 3 and 8. 
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Discussion of 
“SUPPRESSION OF THE FLUID-INDUCED VIBRATION OF CIRCULAR 
CYLINDERS” 


by Peter Price 
(Proc. Paper 1030) 


PETER PRICE.!—The author wishes to thank Mr. Housner for his inform- 
ative discussion of the response of a stack to intermittent forces. Two points 
in particular raised by him prompt further discussion. They are: (1) That 
the stack vibration observations made by him suggest that the flow patterns 
around actual stacks and around the wind tunnel models may be different, and 
(2) the influence of damping on the water channel results quoted in the paper. 

It is the writer’s belief that the observed stack vibration is reconcilable 
with the wind tunnel results if the level of turbulence in the air streams ap- 
proaching the stack and the wind tunnel model are kept in mind. The effect 
of turbulence on the drag coefficient (which reflects the character of flow 
around the cylinder) has been studied extensively.2 It has been found that in 
wind tunnels where the turbulence level is high the critical Reynolds number 
range may be in the region of 60,000 to 200,000,but may be from 200,000 to 
600,000 in a low turbulence tunnel. The 120,000 to 300,000 range reported in 
the paper reflects the high turbulence level of the Stanford tunnel, caused by 
a honeycomb located just upstream of the test region. 

The turbulence level in a wind tunnel stream is bound to be higher than 
will be encountered in normal natural airstreams and thus the critical range 
of Reynolds numbers for a real stack will be higher than wind tunnel tests 
would indicate. This implies that supercritical flow may not be established 
around a stack until the Reynolds number approaches 1,000,000. 

In this light it seems not unlikely that the stack observed by Mr. Housner 
was operating in the critical range (Reynolds numbers of 460,000, 663,000, 
and 1,020,000 for wind speeds of 18, 26 and 40 fps respectively). The larger 
amplitude at the 26 fps wind speed would then be consistent with the wind 
tunnel results where critical range vibration amplitude was larger than that 
encountered in the low supercritical regime. The frequency observations 
made by Mr. Housner for the 18 fps condition also reinforce this belief. It 
would be interesting to test the effect of boundary layer trips or extensive 
surface roughness (devices which artificially produce the turbulent boundary 
layer typical of supercritical flow) on the vibration amplitude of the stack 
mentioned. If the stack has been operating in the critical regime the arti- 
ficially induced supercritical flow conditions should result in a reduced vibra- 
tion amplitude. 


The effect of damping on plain cylinder vibration has been discussed by 


1. Athlone Fellow, Royal Aircraft Establishment, Farnborough, Hampshire, 
England. 


2. Goldstein, S: Modern Developments in Fluid Dynamics, Oxford, The 
Clarendon Press, Vol. HI, 1938. 
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Meier -Windhorst.? Damping reduces vibration amplitude substantially and 
vibration frequency very slightly. If hydrodynamic forces and the fraction of 
critical damping (c) are constant the product of spring stiffness (k) and maxi- 
mum double amplitude (B)max should be constant and the curve in Figure 3 
should be a straight line. The fact that the run 4 (I = 0.335 in.-lb.-sec. 2) 
produces a KBmax factor much lower than that of runs 7 - 11 (I = 2.151 to 
6.722 in.-lb—sec.2) and a maximum amplitude below the projected straight 
line through the Bmax values for runs 7 - 11, would indicate that the damping 
was greater for run 4 than runs 7 - 11, rather than as Mr. Housner suggests. 
The proportion of the beneficial effect of the shroud due to increased 
damping was not directly measured in the water channel tests, which were of 
an exploratory nature, made solely to pick the most likely form of suppres- 
sor. A measure of the importance of damping can be obtained, however, from 
the information plotted in Figure 5. If the very pessimistic assumption is 
made that the entire beneficial effect of the shroud at low k values is due to 
an increase in damping 


“shroud B plain cyl. 


“plain cyl. B vax shroud 


This ratio has been computed throughout the k range and is as follows: 


“shroud 

“plain cyl. 
30 1.84 
80 1.81 
150 1.68 
200 1.58 
300 2.14 
400 5.50 


20 


Since the vibration frequency at maximum amplitude is constant as k is 
varied, the tip velocity of the vibrating cylinder will diminish as maximum 
amplitude decreases, with a consequent reduction in damping. It can be seen 
from the table above that the results predict an increase in damping as k in- 
creases and amplitude decreases. It thus becomes apparent that the dimin- 
ished amplitude is the result of a reduction in the driving force. This re- 
duction is the result of the changed character of the flow around the cylinder 
caused by the shroud. It is significant that the effect is most marked in the 


cases where double amplitudes are 0.7 diameters or less, the region of great- 


est interest. 


3. Meier-Windhorst, A.: Flatterschwingungen von Zylindern im 
gleichmassigem Flussigkeitsstrom. Mitt. des Hydraulischen Institute der 
Technischen Hochschule Munchen, Heft 9, Munich, 1939. 
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Discussion of 
“SOME OBSERVATIONS ON OPEN CHANNEL FLOW AT SMALL 
REYNOLDS NUMBERS” 


by Lorenz G. Straub, Edward Silberman, and Herbert C. Nelson 
(Proc. Paper 1031) 


CORRECTIONS. —The figures to which Prof. Posey refers in his discus- 
sion of this paper in Proc. Paper 1154 appear on page 1154-17. 


LORENZ G. STRAUB,* M. ASCE, EDWARD SILBERMAN, ** and HERBERT 
C. NELSON, *** A.M. ASCE.—The comments of Professors Lansford and 
Robertson and Professor Posey augment the paper and their efforts are ap- 
preciated. Little additional comment is required. 

Regarding the filament of maximum velocity in laminar flow, the measure- 
ments described in the paper do not permit making a positive statement. 

Since abou 15 minutes were required for each pitot tube reading, it was not 
feasible to search for the filament of maximum velocity. Velocities were ob- 
tained by the tube at no closer than 0.2 in. below the surface (about 0.25 depth) 
and at 0.25 in. intervals below this. Within these limitations, no velocity was 
found which exceeded the surface velocity in laminar flow. On the other hand, 
for most turbulent velocity measurements, with the same limitations, a maxi- 
mum was found beneath the surface, even at Reynolds numbers as low as 
6,000. Since the original paper was published, another thesis has been com- 
pleted in which laminar flow velocities were measured by an electrical con- 
ductivity method.(1) The closest readings beneath the surface were about 0.2 
depth and, again, no maximum was found beneath the surface. It appears, 
therefore, that the maximum velocity in laminar flow may be at the surface. 
(The isovels in Professor Posey’s figure from Pickett’s measurements seem 
to indicate maximum velocity at the surface, also.) 

The conclusions regarding the applicability of laminar flow theory to flow 
in open channels were intended to apply mainly to friction factor and to wall 
shear stress—conclusions 1 and 3, page 20. Professor Posey is correct in 
stating that the predicted velocity profiles are not accurately comparable 
with theory (at a distance from the wall). It is believed that the discrepancies 
are surface effects, probably attributable to surface tension. There is every 
reason to believe that the water and kerosene used in the experiments des- 
cribed in the original paper behaved as Newtonian fluids and that the discrep- 
ancies are not explainable by use of differential equations based on non- 
Newtonian fluids. Had Professor Posey included mean flow data with his 
figure taken from Pickett’s measurements, a better comparison could have 
been made with the theoretical results (given graphically in Professor Posey’s 
first figure and analytically in Appendix B-4). 


* “Director, St. Anthony Falls Hydr. Lab., Univ. of Minnesota, Minneapolis, 
Minn. 
** Associate Prof., St. Anthony Falls Hydr. Lab., Univ. of Minnesota, 
Minneapolis, Minn. 
*** Hydr. Engr., Kimberly-Clark Corp., Neenah, Wis. 
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In connection with Fig. 4 of Lansford and Robertson, the following tabula- 
tion of theoretical laminar wall shear distribution for the same channel 
(taken from Appendix B-3) may be interesting: 
LAMINAR WALL SHEAR DISTRIBUTION 
90-Degree Triangular Channel 


/1.414d To/To £/1.414d To/To 
0.05 0.424 0.25 1.135 
0.10 0.700 0.30 1.210 
0.15 0.868 0.40 1.323 


0.20 1.014 
Symmetrical about §/1.414d = 0.5 


0.50 


The following corrections should be made in Appendix B of the paper: 


Page 24 - No. 2 - In the formula for w, the fraction after the summation 
2n + 1 ab 

20 
2n+1mb 


5) qa Also, the first fraction in the equation for K should contain 


sign should contain in the denominator the factor cosh rather than 


rather than a in the denominator. 


Page 25 - No. 3 - In the formula for w, a bracket should enclose the en- 
tire expression following d2, 


2 
Page 25 - No. 5 - The factor (1 + 2 ) should be eliminated from the 
d 


numerator and denominator. 

Page 26 - No. 6 - In the lower figure only, the coordinate axis should be 
rotated 90 degrees counter-clockwise, so that the x-axis is vertical, positive 
upwards, and the y-axis is horizontal, positive to the left. 

Page 26 - No. 7 - The superscript over q in the heading of the table refers 
to the List of References and should be 18 rather than 21. 


Page 27 - No. 8 - The formula for w should be w = 7, b7H. 


Page 28 - No. 9 - The formula for w should be w = aon d7H. 
Page 28 - No. 10 - The reference for w should be 19 rather than 22. 
In Fig. 4 of the discussion by Lansford and Robertson, the ordinate scale 


should contain T, rather than Ty to be consistent with Appendix A. 


REFERENCES 


1. Howe, H. “Laminar Flow Velocity Distribution in a 120° Triangular Chan- 
nel,” Unpublished M. S. Thesis, University of Minnesota, Minneapolis, 
Minn., 1956, 143 pages. 
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Discussion of 
“THE BOUNDARY LAYER DEVELOPMENT IN OPEN CHANNELS” 


by J. W. Delleur 
(Proc. Paper 1138) 


T. BLENCH* M. ASCE.—With a properly designed parallel throated long- 
crested weir equipped with suitable gauging in the channel upstream, the dis- 
charge equation should be deducible to about 0.5% accuracy if the engineer 
will observe the velocity distribution at the gauging section and in the vicinity 
of the critical section in the throat. The equation to use is that expressing 
the fact that, somewhere near the end of the throat, depth must be critical as 
expressed strictly by making the energy of flow (specific depth, specific head, 
still-pond depth) a minimum after allowing for the effect of the kinetic energy 
correction factor (Ref. A). The principal source of error will then be in es- 
timating the friction loss between the gauging and critical sections. Suitable 
designs have been in use extensively in the Punjab Irrigation canals since as 
early as 1927. The present paper interests the writer because it draws atten- 
tion to the need for observing velocity distribution near the critical section 
for every accurate rating of a long-crested weir and suggests that further re- 
search might permit the velocity distribution to be deduced. But with great 
accuracy demanded by the practising engineer the reliability of deduction 
must be fully established. The following queries, with their implied criti- 
cisms, are made with this practical object in view; the writer pretends to no 
expertness in pure boundary layer theory and must expect at least some of his 
questions to appear rather elementary. 


i. Does Fig. 2, or its related theory, indicate that a normal velocity dis- 
tribution could not be established in the throat of a long-crested weir, no mat- 
ter how long the crest might be? If so can the theory be accepted as even 
approximately correct? If not what does a turning value of eta of 0.13 repre- 
sent practically (see Figs. 2, 8)? 


ii. What is H? Is it \ q2/g or y/ka2/g or the value of depth that makes 
energy of flow, viz h + hq2/2g a minimum allowing for k being variable, where 
k is kinetic energy correction factor? 


iii. Does the paper apply to anything else but critical, or very nearly 
critical, conditions? If so, what is the justification for equation (15) which ap- 
pears to compel critical condition? 


iv. Since the phenomenon under investigation is that of the boundary layer, 
that owes its existence to friction loss and is essentially turbulent, why is 


friction loss ignored generally and why are the Reynold’s stresses ignored in 
equation (4)? 


v. Specifically, concerning equation (4), is it merely an abbreviated way of 
writing the straightforward momentum equation allowing for pressure gradient 


* Prof. of Civ. Eng. and Cons. Engr., Univ. of Alberta; Edmonton, Alta., 
Canada. 
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and replacing pressure gradient by its equivalent from the standard hydro- 
dynamic equation of motion which assumes non-pulsating (and therefore non- 
turbulent) flow (Ref. 3 Chapter IV. Sec. 52)? If so, why is this form used 


when there is no pressure gradient but there is turbulence which requires the 
introduction of Reynold’s stresses Pu'v' etc.? 


vi. In equations (18) sigma appears to be H3/2 divided by itself yet it is 
used as a variable. How is this explained? 


vii. Has the virtual channel the same kinetic energy and momentum per 
unit weight of fluid as the actual channel? If not what is the justification for 
using it? 


viii. Why is single log paper used for plotting velocity distributions in 
Figures 5-7 when a seventh root law is used in equation (8)? 


REFERENCE 
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Discussion of 
“WIND INDUCED VIBRATION OF CYLINDRICAL STRUCTURES” 


by Joseph Penzien 
(Proc. Paper 1141) 


T. AU,* A.M. ASCE.—The investigation on vibration of steel stacks has 
been activated by the many instances of reported excessive vibrations of 
welded steel stacks in recent years. A better understanding of the subject is 
still essential in order to design such structures adequately. The paper has 
indeed contributed some additional information much needed in this field. 

Below the critical range of Reynold’s numbers, the vibration of stacks in a 
direction normal to the wind is generally attributed to the von Karman vortex 
trail, and the theory has been validated by experimental results. At or above 
the critical range, hypotheses have been advanced, perhaps without sufficient 
evidence in some instances, to predict the behavior of vibrations. One of the 
erroneous conclusions is that since the shedding of vortices become a- 
periodic beyond the critical range, the turbulent wake would inhibit vibra- 
tions. This obviously contradicts some most recent observations reported in 
literature in which excessive vibrations of stacks involve Reynold’s numbers 
considerably above the critical range. The wind tunnel investigation of the 
dynamic behavior of cylindrical sections under the action of wind at critical 
range reported in the paper further confirms the field observations. 

A recent investigation of vibrations of full-size stacks by Ozker and 
Smith(1) has been cited by the author to compare with his own results. Two 
pertinent points in the former work related to this paper may be elaborated 
here: 


1. The stack vibrates at its own natural frequency at all times regardless 
of wind velocity probably for all conditions above critical range of Reynold’s 
numbers. The alternating force is automatically resonant with the natural 
frequency of vibration. This concept is now accepted in a wider circle al- 
though the exact mechanism behind the formation of the alternating exciting 
forces is still a matter of conjecture. As pointed out by Sherlock, (2) for 
vibration of stacks in wind above critical range, the von Karman theory has 
served a useful framework on which to hang coefficients based on experience, 
but it can not be said that the reality of the von Karman trail has been estab- 
lished. The results from wind tunnel tests of cylinders presented in the paper 
seem to strengthen this impression. 


2. The displacement of the top of stacks increases with increasing wind 
velocity. In other words, once the vibration is initiated, the amplitude of 
vibration will be increased continuously, following approximately a parabolic 
curve, as the wind velocity is increased. This conclusion seems to be at 
variance with the experience of others. Since similar results have been ob- 
served from the stacks in several plants of the Detroit Edison Company as 
reported by Ozker and Smith, there is no reason to believe that it is merely 


* Asst. Prof. of Eng. Mech., Univ. of Detroit, Detroit, Mich. 
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a coincidence. Unfortunately, the paper does not shed any new light on this 
controversy since the range of Reynold’s numbers in the wind tunnel tests is 
between 105 and 106 while that of vibrating stacks may be almost as high as 
107. The pulsating response as indicated by the maximum dynamic reaction 
in the paper nevertheless casts some shadows on the above conclusion. 

The Strouhal number and the drag coefficient for rigid cylinders of 4-inch 
and 12-inch diameters have been found by Delany and Sorensen(3) for 
Reynold’s number R up to 2.3 x 106. Between R = 10° and 106, the results 
are remarkably similar to the S-curve of Relf and Simmons and the Cp-curve 
of Gottingen shown in Fig. 5 of the paper. For R greater than 106, the follow- 
ing values are given: 


R Cp s 
106 0.40 0.38 
1.5x106 0.43 0.42 
2.0x106 0.50 0.35 


Note that the value of Cp is increasing with R up to the limit of the range but 
the curve of S dips for R over 1.5 x 106, 


REFERENCES 
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3. N. K. Delany and N. E. Sorensen, “Low-speed Drag of Cylinders of 
Various Shapes” NACA Technical Note 3038, Nov. 1953. 
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Discussion of 
“A PRESSURE LINE CONCEPT FOR INELASTIC BENDING” 


by Frank Baron 
(Proc. Paper 1157) 


JOHN A. HRIBAR,! J.M. ASCE.—The author has presented here a “tool” 
that can become quite useful in quickly estimating the effects of plasticity on 
the behavior of structural elements subjected to flexural loads. The contri- 
bution hinges on the pictorial design procedure that is used. If this procedure 
is used to obtain a quick result, it results in only an estimate unless one is 
experienced in the method and has a keen sense of judgment in balancing 
statics pictorially. Of course only an estimate may be the desired result. 
But for an accurate result several “checks” would have to be made to obtain 
the best “fit.” 

This procedure may neither be more rapid nor more accurate than 
present methods used on certain materials and cross-sections. Two ques- 
tions are raised. Could not members of steel, that has an easily approxi- 
mated stress-strain curve and comprises the major portion of structural 
material, be analyzed by present methods with the same ease? Could not the 
ultimate load on a structural element of any material and any cross-section 
be as easily and more accurately calculated by the theory of plasticity? The 
writer has raised these questions to emphasize that the theory of plasticity 
may be as rapid, more accurate and less artificial in some cases. 

The writer agrees that as the cross-section or the stress-strain curve 
becomes more unwieldy the “tool” becomes more powerful. Also the by- 
products of the procedure are quite useful in a thorough analysis of a struc- 
ture. The pressure line concept will speed up estimates on structural ele- 
ments of this type. This is a contribution that is well received in a time 
when there is much discussion of ultimate or plastic design. 


Graduate Teaching Asst., Dept. of Civ. Eng., Carnegie Inst. of Technology, 
Pittsburgh, Pa. 
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Paper 1322 


Journal of the 
ENGINEERING MECHANICS DIVISION 


Proceedings of the American Society of Civil Engineers 


DISPERSION OF LONGITUDINAL WAVES* 


by E. Volterra,** M. ASCE 
(Proc. Paper 1322) 


SUMMARY 


The problem of dispersion of longitudinal waves in elastic rods of infinite 
lengths and of rectangular cross-section is discussed by applying the one- 
dimensional theory of wave propagation based on the Method of Internal Con- 
straints. The results given by this theory are compared with those given by 
the Elementary, Love, Bishop approximate theories, and, in the case of two- 


dimensional elasticity, with the exact theories given by Lord Rayleigh and 
Lamb. 


NOMENCLATURE 
The following nomenclature is used in the text: 

2a, 2b = lateral dimensions of the cross-section of the bar 
c = phase velocity of longitudinal waves in a bar 

= group velocity of longitudinal waves in a bar 

E 
plate velocity 
= velocity of Rayleigh surface waves 


5 = velocity of longitudinal waves of infinite wave-length in a 
bar (bar velocity) 
= (is o\(1-20)P = velocity if dilatation waves in an unbounded 
medium 


Note: Discussion open until December 1, 1957. Paper 1322 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 3, July, 1957. , : 

The material in the paper is based on an investigation, conducted at the 


Rensselaer Polytechnic Institute, Troy, N. Y., under the sponsorship of 
the Office of Ordnance Research, Contract No. DA-30-115-ORD-709. Pro- 
ject No. 454.13. 

** Prof. of Eng. Mechanics, The Univ. of Texas, Austin, Tex. 
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“Vopdea) ~ velocity of distortion waves in an unbounded medium 


= Young’s modulus 


oO = Poisson’s ratio for the material of the bar 
Pp = density 


Ty,Tz = radii of gyration of the cross-section of the bar with respect to the 


y and z axes 
=\/r24 r2 


y = polar radius of gyration of the cross-section of the 


bar with respect to the x axis 
time 
= components of the elastic displacement in the x, y and z axes 


= rectangular coordinates 


= wave number 


INTRODUCTION 


The paper discusses the problem of the dispersion of longitudinal waves in 
elastic straight rods of infinite lengths and of rectangular cross-sections by 
applying the one-dimensional approximate theory of wave propagation based 
on the “Method of Internal Constraints,” which theory was presented in pre- 
vious papers. (1,2,3) 

This theory is based on the assumption that during motion the sections 
originally normal to the axis of the bar remain plane. 

Considering a bar of any cross-section shape referring it to a system of 
orthogonal cartesian coordinates x, y, z, where x coincides with the axis of 
the bar and y and z with the principal axes of inertia of the cross-section of 
the bar, the above condition is satisfied by supposing that the components u, 
v, w of the elastic displacement along the x, y and z axes are expressed by 
the following equations: 


v(x,t) 


y A(x,t) 


z u(x,t) 


where v,A andy are three functions a priori unknown of the variable x and 
of the time t. 

On the assumptions that no external forces (body and surface forces) are 
applied to the bar, that the inertia characteristics of the cross-section of the 
bar are constant, that the bar has infinite length, the following equations of 
motion are derived (see Appendix to the paper): 


1322-2 
c 
E 
r 
Pp i 
u, Vv, W 
2 
20a 
° 
a 
7 
(1) 
a. 


(1-20) 3x 


In equations (2) the quantities r, and rz represent the radii of inertia of the 
cross-section of the bar, E,o and P respectively Young’s modulus, Poisson’s 
ratio and the density of the material of the bar. From equations (2) the dis- 
persion curves, i.e. the relationship between the phase and group velocities of 
propagation of longitudinal waves c and cg versus the wave-length L, can be 
derived. 

Dispersion curves are in general obtained by applying the appropriate 
boundary conditions to the general equations of motion of the medium, and 
exact solutions can be derived only in particular cases among which are the 
case of an infinite long rod of circular cross-section and the case of a plate 
of infinite extent. The first problem was investigated in terms of the general 
equations of elasticity by L. Pochhammer in the year 1876.(4) An outline of 
this exact theory is discussed by A. E. H. Love(5) and numerical values of 
dispersion curves have been calculated by D. Bancroft and R. M. Davies.(6,7) 

The second problem was investigated by Lord Rayleigh in the year 1888, (8) 
by H. Lamb in the year 1917,(9) by R. E. D. Bishop in the year 1953.(10) 
Numerical values of dispersion curves have been calculated by J. P. 

Marsden. (11) 

Instead, the problem of wave propagation in an elastic bar of rectangular 
cross-section in terms of the equations of the Mathematical Theory of Elas- 
ticity has not been solved and therefore approximate solutions must be sought. 

The simplest theory of longitudinal waves in elastic bars predicts that 


such waves travel with a constant velocity cy -V ; (bar velocity), and there- 


fore without dispersion. (12) The corresponding equation of motion is: 


_ yg (3) 
p 
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If lateral inertia of the elements of the bar is considered, the terms - 


4 
o2 r? ; Se must be added to equation (3) to take into account this effect, 
P 0,20,2 


rp representing the radius of gyration of the cross-section of the bar about 
the x axis. The differential equation of motion (3) is modified to the form 
given by A. E. H. Love in 1934;(13) 


at? p dx? P 3x7dt? 
If the shear stresses are taken into account the equation of motion takes the 
form given by R. E. D. Bishop in 1953:(14) 


dt? Bx? P 9x7at? P 2(1+0)p dx* 


Bishop(14) has improved the results of equation (5) by the method of succes- 
Sive approximations. In the Bishop improved theory the geometrical form of 
the cross-section of the bar has to be introduced. In the particular case of 

a bar of rectangular cross-section (2a x 2b) the equation of motion becomes: 


If in equation (6) the terms associated with shear stress are neglected, the 
Love improved equation is obtained for a bar of rectangular cross-section as 
derived by Bishop under the following form: 


| 
at? ax? 3 at? 
3 
6(1+c) p ox* 
20° (+0) p [ Wi 
3 
(6) 
2,2 7 vil 
6 
4.44 x6 
6(l+0) p 5 3 dx 
ail 
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In_ the following text the dispersion equations which will be obtained from the 
above equations will be referred to as derived from: 


Elementary Theory from equation 3 
Love Theory from equation 4 


Bishop Theory from equation 5 
Bishop Improved Theory from equation 6 


Love Theory Improved by Bishop from equation 7 
Theory Based on the Method of Internal Constraints from equation 2 


The object of the present paper is to discuss the problem of dispersion of 
longitudinal waves in elastic rods of infinite lengths and of rectangular cross- 
sections by applying the approximate one-dimensional theory of wave propa- 
gation based on the ‘Method of Internal Constraints” and to compare the re- 
sults with those given by the other approximate one-dimensional theories of 
wave propagation and with the exact theories which are deduced from the 
Mathematical Theory of Elasticity. 

The results of such an analysis clearly show that in the two numerical 
cases considered, i.e. the case of a bar of square cross-section and the case 


of a rectangular bar (2a x 2b) having ratio B = 2 = 0.125, the use of the theory 


based on the “Method of Internal Constraints” gives the closest values to the 
exact solutions, which is in accordance with experimental findings. In fact, 
according to recent experiments by R. W. Morse(15,16) on dispersion effects 
in bars of rectangular cross-sections, the dispersion curves of a bar of 
Square cross-section should give very close results to those obtained with a 
bar of circular cross-section if the ratio of the diameter of the cylinder to the 
side of the square cross-section is 1.13. 

The dispersion curves for a bar of rectangular cross-section should be in 
good agreement with the theoretical plane strain case when the ratio of the 


lateral dimensions of the cross-section of the bar B = " < 0.125. 


The Phase-Velocity of Longitudinal Waves 
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where A is a constant, ¥ = a the wave number, L the wave-length, c the 


phase velocity andi=)y-1. By substituting equation (8) into equations (3), 
(4), (5), (6) and (7) one finds the following expressions for the phase velocity: 
In the Elementary Theory: 


Y=%=1 (9) 
c 


In Love Theory: 
3 


(1 +8?) 6" 


In Bishop Theory: 


In Bishop Improved Theory: 
3 


L 2(1+0) 6(1+0) 


122 (14p*) - 8? | 
3 


In Love Improved Theory: 


¥*64r* 2(1+0)F + Y? [ (148?) + 1| 
3 


= A 
ely(x=ct) 


ety (xect) 


re 

0 ao) 

1? o? (1+87 ) 8? 
22 2\p2 
(1+8 )5 
1 + 4 4 
3 
= 
= 
(12) 
0 

where 

Let: 
a4) 
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where A, B and C are three constants. 

Substituting equations (14) into equations (2) one obtains three homogene- 
ous algebraic equations in A, B and C whose determinant set equal to zero 
yields the following frequency equation: 


2 


(27)*876? 


12(5a+6) (1+6?) (27) 767 +144(4a+4) 
(21) *B76* 


+2a+5 | + 


4 [ 12(1+a) (1+87) (27)*6*+144( 3a+2) | 


2(1+0)° (2) 


a+2 
2(1+0)? 


In equations (9), (10), (11), (12), (13), and (15) 


Y 


6 = 2 
L 


From the above equations it follows for the two limiting cases of very 


long (6 > 0), or very short waves (§ > °°), the following valucs for the phase 
velocities: 
In the Love Theory: 


In the Bishop Theory: 


(velocity of distortion 
= Cy waves in an unbounded 
2(1+0)p medium) 


Bate 
— = 
= 
r 
2(1l+0)? 
a 
4 
= 
4 4 ; 
20 
1-260 
= 
3 for Yel 4, 
Y=0 c= 0 
q 
1 
2(1+0) 
ie 
/ 
3 :. 
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In the Bishop Improved Theory: 
for 


for 


for 
for 5 
Instead, for 6>°© equation (15) becomes: 


2(1+0) 2(1+0)* 


The roots of equation (16) are: 


leo 
(1+0) (1-20) 


E(1-o) 


(1+0) (1-20) p waves in an unbounded 
medium) 


(velocity of dilatation 


2(1+0) 


E i (velocity of distortion 
ao. waves in an unbounded 
medium) 


In the case of an infinite plate (plane strain case), the relationship between 
the phase velocity c of sinusoidal waves and their wave length L is given by 
the following transcendental equation (7), (8), (9): 


2 


tanh nb _ 4 (17) 
tanh mb (vy? +n”)? 


where Y = 2 is the wave number 


2 
2040) 2(1+0)p 
In the Love Theory Improved by Bishop: 
c= Cy 
— 
c #0 
xg 
ay 4 
a 
2 3 
q 
‘ 
7 


the velocity of dilatation 
waves in an unbounded 
medium 


the velocity of distortion 


waves in an unbounded 
2p(1+0) medium 


2b the thickness of the plate 


From equation (17) it follows that for very long waves (6 = 2 > 0) 


(plate velocity) 


for very short waves (6>~): c> Cg velocity of Rayleigh 
surface waves (x)(10) 


In the case of a very thin beam (plane stress case), the dispersion equation 
is given by the following transcendental equation (17), (10): 


tanh pa _ 4 y* pq (20) 
tanh qa 


velocity of distortion 
waves in an unbounded 
medium 


the plate velocity 


the depth of the bar 
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From equation (20) it follows that for very long waves 


bar velocity 


for very short waves equation (20) reduces to the cubic equation in ( = 2 


similar to the equation for Rayleigh surface waves; (10) 


22 
+ ( = 16 - 
p 


- 16 (1-4) =0 
“p 


For o = 0.29 the principal root of equation (22) is: 


& = 0.9154 0.5699 


Co fy 


Equations (17) and (20) have been solved numerically by J. P. Marsden(11) of 
the Department of Physics at the University College of Wales in Aberystwyth, 
for values of Poisson’s ratiog = 0.29.* 


In figure 1 the dispersion curves — versus the ratio 2a are given fora 
Oo . 
bar of rectangular cross-section with ratio B = ° = 0.125. The results for the 


lowest mode according to the theory based on the “Method of Internal Con- 
straints” (equation 15) are compared with the results given by the plane 
stress case (equation 18), by the Elementary Theory (equation 9) and by the 
Love and Bishop simple and improved theories (equations 10, 11, 12 and 13). 


In figure 2 the dispersion curves = versus the ratio 2 are given for a 
bar of rectangular cross-section with ratio 8 = = 0.125, according to the 
theory based on the “Method of Internal Constraints” (equation 15). 


Figure 3 gives the dispersion curves = versus the ratio i between 


* The author is indebted to Professor R. M. Davies, Head of the Department 
of Physics at the University College of Wales in Aberystwyth, for having 
furnished him with these numerical data which have been used for compari- 

son in some of the figures enclosed in the text. 
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radius of gyration of the cross-section and wave-length for a bar of square 
cross-section. The curves obtained from equation (15) in the case 8 = 1 are 
compared with the curves calculated by R. M. Davies(7) using the Pochham- 
mer’s solution for a cylindrical bar. . 

In figure 4 the results for the lowest mode of the dispersion curves - 

versus the ratio z for a bar of square cross-section according to the theory 
based on the “Method of Internal Constraints” (equation (15) are compared 
with the results given by the other theories, i.e. with the lowest modes of the 
Pochhammer’s solution according to R. M. Davies, with Love, and Bishop 
simple and improved theories. In the computations of these dispersion 


curves the value 9 = 0.29 has been assumed for Poisson’s ratio of the mate- 
rial of the bar. 


The Group Velocity of Longitudinal Waves 


The velocity with which the energy of a pulse is propagated is called its 


group velocity and is given in terms of the phase velocity c and of the wave- 
length L by the equation: 


(23) 


From equation (23) it follows that the group velocity of longitudinal waves is: 


In the Elementary Theory: (24) 


In Love Theory: (25) 


In Bishop Theory: 


3 6(1+0) 
In Love Theory Improved by Bishop: 


4 ( & )2(140)F + mo (1+p? )6? 


2,42 
4 ( & )2(140)F 4 (14+82)6241 
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In Bishop Improved Theory: ‘ 
df, (5) df,,(6) df.(6) 
cy Coke) a6 ¥ 
cy =clile 
4 ( )* +2 ( & )? £,(6) 
Cc 
where: 
f,(6) = 2(1+0) 7* F 
=- 27* F 6% + (1467)6? +1 
df, (5) 
d6 
(6) 
2,2 
ete) 3 
af. (6 
3 
(1+0) 3(1+0) 
F = isd (1+8*) - 2 | 
3 5 
In the Theory based on the “Method of Internal Constraints:” 
| af, (6) df,(6) df.(6) 
d6 d6 a6 
ec { } (29) 
6 4 ( +2 ( )? 
where: 
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\ 2 
67 2(1+40) 


f5(6) = _3(50+6) 1 + a+1 
4(1+0)?r? 6? (1+0)?1r* B? 4(1+0)? 


af, 148? 2 
d6 (1+0)7? B? 6? 


148% 1 _36(e+1) 1 2 


2(l+0)?7? B? 56? (1+0)27* 64 


- 1 _ 18(30+2) 1 1 


(1+0) B? 62 B= 


20 
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In figure 5 the dispersion curves a versus the ratio = 
of rectangular cross-section with ratio B = ° = 0.125. The results for the 


are given for a bar 


lowest mode according to the theory based on the “Method of Internal Con- 
straints” (equation 29) are compared with the results given for the plane 
stress case (as computed by J. P. Marsden), by the Elementary Theory (ac- 
cording to equation 30) and by the Love and Bishop simple and improved 
theories (according to equations (25), (26), (27) and (28). 
c 
In figure 6 the dispersion curves = versus the ratio 2 are given fora 

bar of rectangular cross-section with ratio B = : = 0.125, according to the 


theory based on the “Method of Internal Constraints” (equation 29). 
In figure 7 the dispersion curves = versus the ratio z (between radius 

of gyration of the cross-section and wave-length) are given for a bar of 
Square cross-section. The curves obtained from equation (29) in the case 
B = 1 are compared with the curves calculated by R. M. Davies(7) using the 
Pochhammer solution for a cylindrical bar. 


c 
In figure 8 the results for the lowest mode of the dispersion curves = 


fe) 
versus the ratio t for a bar of square cross-section according to the theory 


based on the “Method of Internal Constraints” (equation 29) are compared with 
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the results given by the other theories, i.e. with the results for the lowest 
mode of the Pochhammer solution according to R. M. Davies, and of the Love 
and Bishop simple and improved theories (according to equations (25), (26), 
(27) and (28) . 


APPENDIX 


The Equations of Motion for Longitudinal Waves Derived by the 
“Method of Internal Constraints” 


The Strain Energy per unit volume of the bar will be expressed by: 


y = —E(i-¢) 


| 
igs 
| z 
ox 
The Total Potential Energy of the bar will be expressed by: 
u-f vav--f af v dydz (32) 
» 2 A 
where V is the volume of the bar, ¢ its length and A its cross-section area. 
Substituting equations (31) into equations (30) and (32), and performing the rs 
integrations with respect to the variables y and z, calling: rT 
I= dA 
y A 
| I= Y? dA 4 
A (i 
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the principal moments of inertia of the cross-section, one obtains the follow- 
ing expression for U: 


2(1+0) (1-20) 3x 


OE 
(1+0) (1-20) 


The Kinetic Energy per unit volume of the bar will be expressed by: 


at dt | 


where P is the density of the material of the bar. By substituting into equa- 
tions (33) the expressions for the components of the displacements u, v, w 
given by equations (1), performing the integrations over the cross-section of 
the bar, one obtains the Total Kinetic Energy of the bar expressed by: 


On the assumption that no external forces (body and surface forces) are ap- 
plied to the bar, supposing the bar of infinite length, by applying Hamilton’s 
Principle under the form: 


to 
(U+T) dt =0 


to the whole bar, supposing the displacement to be zero at t, and tg, the equa- 
tions of motion (2) are obtained. 
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B= 0.105 


T= 0.29 


2.0 
2a/L 
Figure 2 
Phase velocity of longitudinal waves in a bar of rectangular 
cross—section according to the Theory of Internal Constraints 
( 3 Modes). 
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0.6 
r/L 
Figure 3 
Phase velocity, Cpe of longitudinal waves of wave length L 
in square cross-section bars of radius of gyration r. 


1, 2, 3, Exact theory (According te R. M. Davies), 
(for cylindrical bars), 3 Modes. 


1', 2', 3", Method of Internal Constraints, 3 Modes, 
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Figure 6 


Group velocity of longitudinal waves in a bar of rectangular 
cross-section according to the Theory of Internal Constraints, 
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Figure 7 


Group velocity, cg, of longitudinal waves of wave length L, 
in square cross-section bars of redius of gyration r. 


1,2 Exact theory (According to R.M. Davies). 1', 2' Theory 
of Internal Constraints. 
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